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ABSTRACT 

We present detailed 16-GHz interferometric observations using the Arcminute Microkelvin 
Imager (AMI) of 19 clusters with L x > 7 x 10 37 W (hso = 1) selected from the Local Cluster 
Substructure Survey (LoCuSS; 0.142 < z < 0.295) and of Abell 1758b, which is in the field 
of view of Abell 1758a. We detect and resolve Sunyaev-Zel'dovich (SZ) signals towards 17 
clusters, with peak surface brightnesses between 5 and 23 <x. We use a fast, Bayesian cluster 
analysis to obtain cluster parameter estimates in the presence of radio point sources, receiver 
noise and primordial CMB anisotropy. We fit isothermal /3-models to our data and assume the 
clusters are virialized (with all the kinetic energy in gas internal energy). Our gas temperature, 
?ami, is derived from AMI SZ data and not from X-ray spectroscopy. Cluster parameters inter- 
nal to r5oo are derived under the assumption of hydrostatic equilibrium. We find the following, 
(i) Different gNFW parameterizations yield significantly different parameter degeneracies, (ii) 
For hjo = 1, we find the classical virial radius, r2oo, to be typically 1.6+0.1 Mpc and the to- 
tal mass Mx(r2oo) typically to be 2.0-2.5X Mx^oo). (iii) Where we have found Mx^oo) and 
^t(^20o) X-ray and weak-lensing values in the literature, there is good agreement between 
weak-lensing and AMI estimates (with M t ,ami/^t,wl = 1 -2^q| and = 1.0 + 0.1 for rsoo and 
r2oo, respectively). In comparison, most Suzaku/Chandra estimates are higher than for AMI 
(with M t ,x/^t,ami = 1.7+0.2 within rsoo), particularly for the stronger mergers, (iv) Compar- 
ison of Tami to T x sheds light on high X-ray masses: even at large radius, T x can substantially 
exceed Tami in mergers. The use of these higher T x values will give higher X-ray masses. We 
stress that large-radius Tami and T x data are scarce and must be increased, (v) Despite the 
paucity of data, there is an indication of a relation between merger activity and SZ ellipticity. 
(vi) At small radius (but away from any cooling flow) the SZ signal (and Tami) is less sen- 
sitive to ICM disturbance than the X-ray signal (and T x ) and, even at high radius, mergers 
affect n 2 -weighted X-ray data more than n- weighted SZ, implying that significant shocking or 
clumping or both occur in even the outer parts of mergers. 

Key words: cosmology: observations - cosmic microwave background - galaxies: clusters - 
Sunyaev-Zel'dovich X-ray - galaxies: clusters: individual (Abell 586, Abell 611, Abell 621, 
Abell 773, Abell 781, Abell 990, Abell 1413, Abell 1423, Abell 1704, Abell 1758a, 
Abell 1758b, Abell 2009, Abell 2111, Abell 2146, Abell 2218, Abell 2409, RXJ0142+2131, 
RXJ1720.1+2638, Zw0857.9+2107, Zwl454.8+2233) 
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1 INTRODUCTION 

The virtues of galaxy clusters are often extolled as, for example, be- 
ing the largest gravitationally bound systems in the Universe, or be- 
ing excellent samplers of the matter field on large scales, or simply 
as being o f fundamental importance to astrophysics an d c osmology 
(see e .g., lwhite et al.lll993LlEke. Cole. & Frenklll996l and lJov et all 
120011) . To make full use of these virtues one needs observations 
that, amongst others things, reach large distances away from cluster 
centres. It would often be very useful to reach the classical virial ra- 
dius « r 2 oo of a cluster, internal to which the average density is 200 
times the closure density. Studying clusters on these scales is im- 
portant for many reasons. First, these measurements are needed to 
characterize the entire cluster volume. Second, they can be key for 
any attem pt at precisio n cosmology, including calibrating scaling 
relations dKaisedll986h . as they are believed to be less susceptible 
to the complicated physics of the core region from e.g., star forma- 
tion, energy feedback from active galactic nuclei and gas cooling. 
Third, the virial radius marks the transition between the accreting 
matter and the gravitationally-bound, virialized gas of clusters and 
thus contains information on the current processes responsible for 
large-scale structure formation. However, there are few such obser- 
vations due to the difficulties of obtaining a signal far away from 
the cluster centre. We now comment on four method s of estimat- 
ing cluster masses (see lAllen. Evrard. & Mantzll201ll for a recent, 
overall review): 

• Spectroscopic measurements of the velocity dispersion of 
cluster members require very high sensitivity at moderate to high 
redshift, and confusion becomes worse as redshift increases and 
as distance on the sky from the cluster centre increas es. Cluster 
masse s ha ve recently been o btained this way in e.g., iRines et all 
d201Qh and lSifon et all d2012h . 

• X-ray observations of the Bremsstrahlung (free-free radi- 
ation) from the intracluster plasma (by convention referred to 
as 'gas') have delivered a great deal of inf ormation on cluster 
physics on a large num ber of clusters (see e.g., Ebe ling et"aT] [ 1998, 
iBonamente et al 120061 and lKotov & VikhlinirJI 20051) . Observations 
are, of course, difficult at high redshift due to cosmic dimming, and 
because the X-ray signal is oc / n 2 f(T)dl, where n is the electron 
density, T is electron temperature, f(T) is a weak function of T, 
and / is the line of sight through the cluster, there is significant 
bias to gas concentration, which makes reaching a high radius dif- 
ficult - however, at low to intermediate redshift there is a small but 
growing number of observations that approach or reach r 2 oo mainly 
with the Suzaku satellite, though the sky background sub traction is 
challenging (e.g jGeorge et aU2009llHoshino et al.ll2010h . 

• Gravitational lensing of background galaxies gives the dis- 
tribution of all the matter in the cluster directly, without rely- 
ing on assumptions obout the dynamical state of the cluster. Any 
mass concentrations along the line-of-sight not associated with 
the cluster will lead to an overestimate of the weak lensing clus- 
ter mass. But the 'shear' signal is proportional to the rate of 
change with radius of the gravitational potential, which changes 
increasingly slowly with radius at large radius, so reaching large 
radius is difficult. Confusion also bears strongly on this diffi- 
culty, and measurement is of course harder as redshift increases. 
Examp le w eak-lensing cluster stud i es inc lude lOkabe & Umetsul 
(120081) and ICorless. King. & Clowd <2009l). for analy s es of indi- 
vidual high-mass clus t ers, an d iMandelbaum & Seliakl d2007l) and 
IRozo. Wu & Schmidtl J201ll) for analyses of stacked lensing pro- 
files for many low-mass clusters. 

• The Sunyaev Zel'dovich (SZ; ISunvaev &Zel'dovich| [l970: 



see e.g lBirkinshawl[T999l and Carlstrom, Holder & Reese 2002 for 
reviews) signal from inverse Compton scattering of the CMB by 
the cluster gas has relatively little bias to gas concentration since 
it is oc J nTdl, and has remarkably little sensitivity to redshift 
over moderate to high redshift; both of these properties make the 
SZ effect extremely attractive. The problem with SZ is that it is 
intrinsically very faint. The first gen eration of SZ telescope s, in- 
cluding the OVRO 4 0-m (see e.g. , iBirkinshaw et all fl984). the 
OVRO 5- m (see e.g.. iHerbig et allll995l) . the OVRO/BIMA ar- 
rays (e.g.. ICarlstrom et aL 19921) and the Ryle Telescope (see 
e.g., iGrainge et al.l 1 19931) had to integrate for a very long time 
to get a significant SZ detection of a single known cluster. The 
new generation, inc luding ACT (s ee e.g.,|Hm cks et al. 201( j and 
Marr iage et al.l201 1[). AMI (see e.g.jAMI Consortium: Zwart et al.l 
12008 ] and | AMI Consortium: Shimwell et alJbOloL AMiBA (see 
e.g jLo. Martin. & ChiuehlbOOll andlWu et alJl2008h. MU STANG 
(see e.g. 
(se e e.g., 



Korngut et al. 201 1 and Mroczkowski et al. 201 1), OCRA 
Browne et alfeOOd and Lancaster et al.l201lL Planc k (see 



e.g iTauber et al.ll2010l and IPlanck Collaboration et aPbOl ll) . SPT 
(see e.g Carlstrom et al. 201 1 and Williamson et al. 201 1) and SZA 
(see e.g., ICarlstrom et al.|[l992l and iMuchovei etalJbOllb are all 
much more sensitive. 

The new generation of SZ facilities include two types of instru- 
ment: ACT, Planck and SPT are instruments with wide fields of 
view (FoV) optimized for detecting CMB imprints in large sky ar- 
eas in a short amount of time - this is a very important ability but, 
for the imaging of a particular cluster, a wide FoV is of no ben- 
efit; in contrast, AMI, AMiBA, MUSTANG, OCRA and SZA are 
designed to go deep and to measure the masses of the majority of 
clusters. 



In lAMI Co nsortium: Zwa rt et al we reported initial 

SZ observations of seven X-ray clusters (selected to have low ra- 
dio flux-densities to limit confusion) that approach or reach r 2 oo- In 
this paper we report on resolved, interferometric SZ observations 
with arcminute resolution that approach or reach r 2 oo in a substan- 
tial sample of X-ray clusters selected above an X-ray flux-density 
limit (plus a radio flux-density limit) and over a limited redshift- 
range (which limits the effects of cosmic evolution); as far as we 
are aware, this is the first time such SZ observations of a large clus- 
ter sample have been undertaken. These measurements are timely 
since complementary lar ge- r X-ray data have recently been ob- 
tained with Suzaku (e.g. jBautz et al.l2009l . lHoshino et al.l2010l and 
iKawaharada et al .1 120 1 Oh . These early Suzaku measurements, de- 
spite the large model uncertainties, are already showing that ICM 
profiles on these scales appear to disagre e with predictions from 
hydrodynamical cluster simulations (e.g., iGeorge et al.ll2009h and 
have drawn attentio n to possible causes such as ICM clumping 
dNagai&Laul201ll) and the breakdown of assumptions such as hy- 
drostatic equilibrium (e.g., Evrardl ll997l) . which can bias the X-ray 
masse s (e .gjRasia. Tormen. & Moscard ini 2004, Meneghetti et al.l 
120101 and lFabian et al.ll201ll) . We stress that SZ observations, like 
those in the optical/IR and X-ray, also have their contaminants and 
systematics, and all four methods are also hampered by projection 
effects. 

Studying large samples of clusters using multiple tech- 
niques is important for building a thorough understanding of clus- 
ter physics. Well-calibrated mass-observable relations are cru- 
cial for current and f u ture cosmological studies - see e.g., 
lAllen. Evrard. & Mantzl d201ll) . To our knowledge, this is the 
largest cluster-by-cluster study for which masses have been derived 
from SZ targeted observations out to the virial radius. The results 
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from this work will be very valuable for detailed comparisons of 
cluster mass estimates. This paper is organized as follows. In Sec. 
2 we describe the sample selection. The data and instrument are in- 
troduced in Sec. 3, while Sec. 4 focuses on the methods applied for 
mapping the data, identifying radio source foregrounds and remov- 
ing them from the maps. The analysis of the cluster + radio sources 
environment is outlined in Sec. 5. Given the difficulty of comparing 
cluster mass measurements from different data, we provide consid- 
erable detail in our results section, Sec. 6. In particular, we present: 
maps; details on the radio source environment towards the clus- 
ters; full cluster parameter posterior distributions internal to two 
overdensities, r 50 o and r 2 oo; an investigation of contaminating radio 
sources (our main source of systematic error); and we compare our 
/?-model parameterization with several generalized Navarro -Frenk- 
White (gNFW) parameterizations. In Sec. 7 we illustrate the ability 
of our methodology to recover the cluster mass even for a clus- 
ter with a challenging source environment. In Sec. 7 we discuss 
our results, in particular, the morphology and dynamical state of 
the clusters and the comparison of SZ-, weak lensing- and X-ray- 
derived cluster masses and large r X-ray and SZ temperatures. The 
conclusions of our study are summarized in Sec. 9. 

Throughout, we assume a concordance ACDM cosmology 
with Q m ,o = 0.3, Q A , = 0.7, Q k = 0, Q bfi = 0.041, w = -1, 
w a = and <x 8 = 0.8. For the probability distribution plots and 
the tables, we take h = H /\00 kins" 1 Mpc -1 ; elsewhere we take 
H = 70 kms -1 Mpc -1 as the default value and also refer when 
necessary to h x = Ho/Xkms' 1 Mpc -1 . All coordinates are at epoch 
J2000. 



2 THE LOCUSS CATALOGUE AND OUR SUB-SAMPLE 

LoCuSS dSmith et all l200l l2004h is a multi- wavelength survey 
of 164 X-ray luminous (L x > 2 x 10 37 W ove r th e 0. 1-2.4 keV 
band in the cluster rest frame febeling et al.l 19981 and lEbeling et al.l 
2000, h 50 = 1) galaxy clusters. The narrow range of redshifts z 
(0.142 < z < 0.295 minimises cosmic evolution. The clusters 
have b een selected from the ROSAT All-Sky Survey dEbeling et al.l 
1 19981. l2000l : iBohringer et all [2004) without taking into account 
their structures or dynamic al s tates. Relevant LoC uSS papers in- 
clude jM^QnTet^(|2009|) and lZhang et~aDd2010h . 

In this work, we study a sub-sample of 19 clusters from the 
LoCuSS catalogue and Abell 175810 (Tab.[Q using 16-GHz inter- 
ferometric AMI data with arcminute resolution. Our sub- sample in- 
cludes only those clusters with 5 > 20° . AMI can observe down to 
lower declinations but suffers from poorer wv-coverage and satellite 
interference at 8 < 20°. We also applied an X-ray luminosity cut, 
L x > 7x 10 37 W (0.1-2.4 keV restframe, h 50 = 1), lower-luminosity 
clusters tend to be fainter in SZ. Contamination from radio sources 
at 16 GHz can significantly affect our SZ detections. For this rea- 
son, we have chosen to exclude clusters with sources brighter than 
lOmJy beam -1 within 10' of the cluster X-ray centre. 

Several studies of the LoCuSS sample of clusters are ongoing. 
These include both ground-based (Gemini, Keck, MMT, NOAO, 
Palomar, Subaru, SZA, UKIRT and VLT) as well as space-based 
(Chandra, HST, GALEX, XMM-Newton and Spitzer) facilities. 
Our AMI SZ data are complementary to other data taken towards 



Abell 1758b was serendipitously observed in the field of view of 
Abell 1758a, a LoCuSS cluster. 



these clusters as they probe the large-scale gas structure, are sensi- 
tive to gas from destroyed density peaks and are particularly ben- 
eficial for obtaining robust cluster masses since the SZ signal has 
long been recognised as a good mass proxy (see e.g.. iMotl et al.l 
l2005h . 



3 INSTRUMENT AND OBSERVATIONS 

AMI consists of two aperture- synthesis interferometric arrays lo- 
cated near Cambridge. The Small Array (SA) is optimized for 
SZ imaging while the Large Array (LA) is used to observe ra- 
dio sources that contaminate the SZ effect in the SA observations. 
AMI's wv-coverage is well-filled all the way down to « 180/1, cor- 
responding to a maximum angular scale of « 10'. AMI is described 
in detail in AMI Consortium: Zwart et al. 

SA pointed observations of all the clusters were taken between 
2007 and 2010 while LA raster observations, which were mostly 
61 + 19 pt hexagonal raster^ centred on the cluster X-ray position, 
were made between 2008 and 2010. Typically, each cluster was 
observed for 20-80 hours with the SA and for 10-25 hours with 
the LA. The thermal noise levels for the SA (cr SA ) and for the LA 
(cr LA ) were obtained by applying the aip^I task imean on a section 
of the map far down the primary beam and free from any significant 
contamination. In Tab. [2] we provide central thermal noise estimates 
for the SA and LA observations; they reflect the amounts of data 
remaining after flagging. A series of algorithms has been developed 
to remove (or 'flag-out') bad data points arising from interference, 
shadowing, hardware and other errors. This is a stringent process 
that typically results in « 30-50% of data being discarded before 
the analysis. A primary-beam correction factor has been applied, 
as the thermal noise level is dependant on the distance from the 
pointing centre. 

The raw visibility files were put through our lo- 
cal data reduction pipeline, reduce , described in detail in 
lAMI Consortium: Davies et all J201(th and exported in fits 
format. Bi-daily observations of 3C286 and 3C84 were used for 
flux calibration while interle aved calibrators sele c ted from the 
Jodrell Bank VLA Survey dPatnaik et all Il992l : iBrowne et aD 
1998; Wilkinson et "aH ll998l) were observed every hour for phase 
calibration. 



4 MAPPING AND SOURCE DETECTION AND 
SUBTRACTION 

Our LA map-making and source-fin ding procedures follow 
AMI Consortium: Shimwe ll et al We applied standard aips 

tasks to image the continuum and individual-channel uvfits data 
output from reduce. 

At 16 GHz, the dominant contaminants to the SZ decrements 
are radio sources. In order to recover the SZ signal, the contribution 
of these radio sources to the data need to be removed; this is done 
as follows. 

2 The observing frequency range given in AMI Consortium: Zwart et al. 
has been altered slightly, as described in AMI Consortium: Franzen et al. 
2010. 

3 A 61 + 19 point raster observation consists of 61 pointings with separa- 
tions of 4', of which the central 19 pointings have lower noise levels, see 
e.g jAMI Consortium: Hurley- Walker et al.N201ll) for example LA maps. 

4 http://www.aips.nrao.edu 
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Table 1. Cluster details. It should be noted that Abell 1758b is not part of LoCuSS. 



Cluster 


Right ascension 
(J2000) 


Declination 
(J2000) 


Redshift 


X-ray luminosity 

/10 37 W 
(/i50 = 1 ; see text) 


Alternative cluster names 


Abell 586 


07 32 12 


+31 37 30 


0.171 


11.1 




Abell 611 


08 00 56 


+36 03 40 


0.288 


13.6 




Abell 621 


08 11 09 


+70 02 45 


0.223 


7.8 




Abell 773 


09 17 54 


+51 42 58 


0.217 


13.1 


RXJ0917.8+5143 


Abell 781 


09 20 25 


+30 31 32 


0.298 


17.2 




Abell 990 


10 23 39 


+49 08 13 


0.144 


7.7 




Abell 1413 


11 55 18 


+23 24 29 


0.143 


13.3 




Abell 1423 


11 57 18 


+33 36 47 


0.213 


10.0 


RXJ1 157.3+3336 


Abell 1704 


13 14 18 


+64 33 27 


0.216 


7.8 




Abell 1758a 


13 32 45 


i CA 11 

+50 32 31 


0.280 


11.7 




Abell 1758b 


13 32 29 


+50 24 42 


0.280 


7.3 




Abell 2009 


15 00 21 


+21 22 04 


0.153 


9.1 




Abell 2111 


15 39 40 


+34 26 00 


0.229 


10.9 




Abell 2146 


15 55 58 


+66 21 09 


0.234 


9.0 




Abell 2218 


16 35 45 


+66 13 07 


0.171 


9.3 




Abell 2409 


22 00 57 


+20 57 50 


0.147 


8.1 




RXJ0142+2131 


01 42 03 


+21 31 40 


0.280 


9.9 




RXJ1720. 1+2638 


17 20 10 


+26 37 31 


0.164 


16.1 




Zw0857.9+2107 


09 00 39 


+20 55 17 


0.235 


10.8 


Z2089 


Zwl454.8+2233 


14 57 15 


+22 20 34 


0.258 


13.2 


Z7160 



• First, the CLEANed LA continuum mapjfl were put through 
the AMI-developed source-extracti on software sourcefind 
(AMI Consortium: Franzen et al to identify and character- 
ize radio sources on the LA maps above a certain signal-to-noise. 
sourcefind provides estimates for the right ascension x s , decli- 
nation y s , flux densit}|E So, and spectral index cfl at the central 
frequency vo for identified radio sources. We impose a detection 
threshold such that we select only those radio sources with a 
flux density ^ 4cr LAp on the CLEANed LA continuum maps, where 
c"la p refers to pixel values on the LA noise maps. The number of 
>4cr LAp sources detected in our LA observations of each cluster is 
given in Tab. O 

• Second, prior to any source subtraction, we run our cluster- 
analysis software, which fits for the position, flux and spectral index 
of the souRCEFiND-detected radio sources using the source parame- 
ters obtained by sourcefind as priors. For some of the less contami- 
nating radio sources, our cluster-analysis software uses delta-priors 
for the source parameters centred at the LA estimates (see Sec. [5] 
for further details) . 

• Third, the source parameters given by the cluster analysis 
were used to perform source subtraction on the SA maps. This was 
done using in-house software, muesli, which is an adaptation of the 
standard aips task uvsub optimized for processing AMI data. The 
flux-density contributions from detected radio sources were sub- 
tracted from each SA channel uvfits file using either the mean val- 
ues for their position, spectral index and flux density derived from 
our Bayesian analysis, when these parameters are not given delta- 
function priors, or, otherwise, using the LA estimates for these 
source parameters. Details of the priors assigned to each of the 



5 The LA continuum maps were CLEANed down to 3cr LA with no boxes. 

6 We catalogue the peak flux of the source, unless the source is ex- 
tended, in which case we integrate the source surface brightness over 
its projected solid angle to give its integrated flux density (see e.g., 
AMI Consortium : Franzen et al 1201 Ol) . 

7 We adopt the convention S oc y~ a . 



sources labelled on the SA maps can be found in Sec. 15.2.21 and 
Tab.H 

• Fourth, after source subtraction, the SA maps were CLEANed 
with a tight box around the SZ signal. In contrast, the LA maps 
and SA maps before source subtraction were CLEANed with a single 
box comprising the entire map. Both the SA and the LA maps were 
CLEANed down to 3o~. 



5 ANALYSIS 

We use our own Bayesian analysis package, McAdam, to estimate 
cluster parameters internal to r 50 o and r 2 oo from AMI data in the 
presence of radio point sources, receiver noise and primordial CMB 
anisotropy. The cluster and radio sources are parameterized in our 
analysis (see below) while the remaining components are included 
in a generalized noise covariance matrix; we note that these are 
the only significant noise contributions because large-scale emis- 
sion from e.g., foreground galactic emission is resolved out by our 
interferometric observations. 

McAdam wa s originally developed by iMarshall et all d2003l) 
and iFeroz et al . (2009) adapted it to work on AMI data . 
The latest McAdam uses MultiNest (Feroz & Hobsonl I2008L 
IFeroz. Hobson & Bri dges 2009) as its inference engine to allow 
Bayesian evidence and posterior distributions to be calculated ef- 
ficiently, even for posterior distributions with large (curved) degen- 
eracies and/or mutiple peaks. This addition has been key to our 
analysis since the posteriors of AMI data often have challenging 
dimensionalities, > 30, primarily as a result of the presence of a 
large number of radio sources in the AMI observations. 



© 0000 RAS, MNRAS 000, 000-000 



SZ study of 19 LoCuSS galaxy clusters with AMI out to the virial radius 5 

Table 2. Observational details. SA and LA noise levels, osa and cr LA > 4<x and the number of sources detected above 4cr LAp on the LA rasters for each 
cluster. Abell 1758b is not part of LoCuSS. 



Cluster o~ sa &LA Number of LA 4<tla sources 
(mJy) (mJy) 



Abell 586 


0.17 


0.09 


23 


Abell 611 


0.11 


0.07 


23 


Abell 621 


0.11 


0.09 


13 


Abell 773 


0.13 


0.09 


9 


Abell 781 


0.12 


0.07 


24 


Abell 990 


0.10 


0.08 


20 


Abell 1413 


0.13 


0.09 


17 


Abell 1423 


0.08 


0.07 


31 


Abell 1704 


0.09 


0.06 


13 


Abell 1758a 


0.12 


0.08 


14 


Abell 1758b 


0.13 


0.08 


14 


Abell 2009 


0.11 


0.14 


18 


Abell 2111 


0.09 


0.07 


22 


Abell 2146 


0.15 


0.06 


15 


Abell 2218 


0.07 


0.10 


15 


Abell 2409 


0.14 


0.05 


15 


RXJ0142+2131 


0.11 


0.06 


22 


RXJ1720. 1+2638 


0.08 


0.10 


17 


Zw0857.9+2107 


0.13 


0.12 


13 


Zwl454.8+2233 


0.10 


0.10 


16 



5.1 Model 



5.2 Priors 



We have modelled the clu ster density profile assuming sph erical 
symmetry using a/?-model (Cavalier e & Fusco-Ferm iano 1978): 



PgO) = 



Pg(0) 



(1) 



where gas mass density p g (r) = jin(r), \i - 1.14m p is the gas mass 
per electron and m v is the proton mass. The core radius r c gives the 
density profile a flat top at low j- and p g has a logarithmic slope of 
3j3 at large ~ . 

We choose to model the gas as isothermal, using the virial 
mass-temperature relation and assuming that all kinetic energy is 
in gas internal energy: 



k B T(r 2 oo) = 



GyuM T (r2oo) 
2r 20 o 

3 V /3 

\4^(200p crit )/ 



M 2 '\r 2m ) 



47r(200p crit ) ; 

M T (r 20 o) 



■ 8.2keV 



\0 15 h- l M Q j \ H J 



2/3 (H(z)\ 2/3 



(2) 
(3) 

(4) 



5. 2. 1 Cluster priors 

The cluster model parameters C = (x c , y c , M T (r 20 o), /g(?20o)>/^ r c , 
z) have priors that are assumed to be separable. x c , y c are the clus- 
ter position (RA and Dec, respectively) and f g is the gas fraction, 
which is defined as 



(5) 



Further details on these priors are given in Tab. 

This set of sampling par ameters has proved sufficient for 
our c luster detection algorithm dAMI Consortium: Shimwell et al.1 
l2010l) and to describe the physical cluster parameters. We empha- 
size that this way of analysing the data is different from the way 
used traditionally, in which an X-ray spectroscopic temperature is 
used as an input parameter. The difficulty with this use of an X- 
ray temperature is that, in practice, the temperature measurement 
usually applies to gas relatively close to the cluster centre (but any 
cooling flow is excised). By sampling from M T and using the M-T 
relation (Eq.[4]), the temperature of each cluster is derived from SZ 
data only and is averaged over the angular scale of the SZ obser- 
vation, which is typically larger than the angular scale of the X-ray 
temperature measurement. This way, although our analysis does not 
yield T(r), it gives and uses a temperature which is representative 
of the cluster volume we are investigating. 



Mj(r 2 oo) and r(r 2 oo) refer to th e total mass and gas temperature 
within r 2 oo (see e.g., Voit 20o3). This relation allows cluster pa- 
rameters within r 20 o to be inferred without assuming hydrostatic 
equilibrium; note that, in our methodology, parameters describ- 
ing the cluster at smaller r (e.g., r 50 o) do, however, assume hy- 
drostatic ea^nlibriuj^ be 
found in AMI Consortium: Rodriguez-Gonzalvez et al. I d2010h and 
lAMI Consortium: Olamaie et al.l j201lh . The good agreement be- 
tween mass estimates from weak-lensing and AMI d ata on 6 clus- 
ters in lAMI Consortium: Hurley-Walker et al.l (1201 ll) supports the 
use of this M-T relation in our analysis. 



5.2.2 Source priors 

Radio sources detected on the LA maps using sourcefind are mod- 
elled by four source parameters, S = (x s , y s , S , a). Priors on these 
parameters are based on LA measurements, discussed in Sec. |4] 

Sources on the source- subtracted SA maps are labelled ac- 
cording to Tab.|4] Delta-function priors on all the source parameters 
tend to be given to those sources whose flux density is < 4cr 5A and 
to those outside the 10% radius of the SA power primary beam. The 
remaining sources are usually assigned a delta-function prior on po- 
sition and Gaussian priors on a and S - However, in a few cases we 
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Table 3. Summary of the priors for the sampling parameters in each model. The value for the redshift and position priors have not been included in this table 
since they are cluster specific. Instead, they are given in Tab.[T]for each cluster. 



Parameter 



Prior Type 



Values 



Origin 



Xc,y c 
P 

M T (r2oo)/^ _1 M 
r c //T ! kpc 

z 

f g (r200)/h- 1 



Gaussian at xx- r 



:60 



uniform 
uniform in log 
uniform 
delta 

Gaussian, cr = 0.0216 



cluster position 
0.3-2.5 
1 x 10 13 - 5 - 5x 10 15 
10 - 1000 
cluster redshift 
0.0864 



Ebeling et al. (2000) 
Marshall et al. (2003) 
physically reasonable, e.g., Vikhlinin et al. 2006 
physically reasonable e.g., Vikhlinin et al. 2006 
Ebeling et al. (2000) 
Larson et al. (2011); Zhang et al. (2010) 



Table 4. Priors on position, spectral index and flux density given to detected 
sources. The symbols correspond to the labels in the SA source-subtracted 
maps. The Gaussian priors are centred on the LA measurements, cr values 
for the Gaussian priors are assigned as follows: for the Gaussian prior on the 
flux-density of each radio source, cr is set to 40% of the source flux density; 
for the spectral index alpha, cr is set to the LA error on a and for the source 
position, cr is set to 60" . 



Symbol 


ncso) 


11(a) 


n(x 5 ,v 5 ) 


+ 


delta 


delta 


delta 


X 


Gaussian 


Gaussian 


delta 


A 


Gaussian 


Gaussian 


Gaussian 



replace delta-function priors on the source parameters with Gaus- 
sian priors as this can increase the accuracy of the source subtrac- 
tion. These wider priors can be necessary to account for discrepan- 
cies between the LA and SA measurements. Reasons for these dif- 
ferences include: a poor fit of our Gaussian model for the power pri- 
mary beam far from the pointing centre, correlator artifacts, source 
variability and source extension. 



SA maps are the FWHMs of the synthesized beams. A 0.6-kA ta- 
per was applied to the SA source- subtracted maps to downweight 
long-baseline visibilities with the purpose of increasing the signal- 
to-noise of the large-scale structure; this typically leads to a « 20% 
increase in the noise. The X-ray images are obtained from archive 
ROSAT and Chandra data. 

We remind readers that when looking at a radio map - neces- 
sarily with a particular uv- weighting - a near-circular image does 
not mean that the SA failed to resolve the SZ signal. Investigating 
angular structure / size requires assessment in wv-space, which can 
be done with a selection of maps made over different uv ranges but 
is optimally done here in uv- space with McAdam. In fact, all the SZ 
decrements in this paper are resolved. 



6 RESULTS AND COMMENTARY 

Out of the 20 clusters listed in Tab. Q] we detect SZ decrements to- 
wards 17. For these clusters we present SA maps before and after 
source subtraction as well as marginalized posterior distributions 
for some cluster and source parameters and mean values of se- 
lected cluster parameters (Tab.|6), with the exception of Abell 2409, 
which was found to have a local environment which renders it un- 
suitable for robust parameter estimation (see Sec. 16.16) . For the 
posterior distributions all ordinates and abscissae in these plots are 
linear, the j-axis for the 1-d marginals is the probability density and 
h is short for /i 100 . It is important to note that, while the posterior 
probability distributions for large-scale cluster parameters reflect 
the uncertainty in the McADAM-derived flux-density estimates, the 
radio source- subtracted maps do not, as they simply use a single 
value (the mean) for each source parameter. The effect of our pri- 
ors on the results has been thoroughly te sted in a previous study by 
AMI Consortium: Olamaie et al.l 120111) . which found that the pri- 
ors used in this parameterization do not to lead to any strong biases 
in the cluster parameter estimates. 

The SA maps have labels indicating the position of detected 
radio sources and their priors (Tab. @]; the square box in these plots 
indicates the best-fit cluster position determined by McAdam. No 
primary-beam correction has been applied to the SA maps pre- 
sented in this paper, unless stated otherwise. The contour levels 
on the SA maps, unless otherwise stated, start at 2<t S a and increase 
linearly from 2 to l0cr SA . On radio-only images, positive contours 
are shown as solid lines and negative contours as dashed lines, but 
on radio+X-ray images, negative radio contours are shown as solid 
lines and X-ray shown as greyscale. The bottom-left ellipses on the 
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Table 5: Source properties for detected sources within 5' of the SZ mean central 
position. The number next to each cluster name denotes the source number; this 
label is used in the plots showing the marginalized posterior distributions for the 
source fluxes. SoMcA is the McADAM-derived best-fit source flux at 16 GHz. a 
is the source spectral index estimated by McAdam and centred at the McADAM- 
derived mean frequency and the last column contains the distance between the 
cluster SZ centroid (as determined by McAdam) and the source. 



Name Right ascension Declination SqMcA a Distance from SZ centroid 





(hh:mm:ss , J2000) 


(°:':" , J2000) 


mJy 






A586_0 


07 


32:20.5 


+31:38:02.8 


0.26 


1.20 


28 


A586_l 


07 


32:19.1 


+31:40:25.6 


0.86 


0.28 


171 


A586.2 


07 


32:11.0 


+31:39:47.6 


0.86 


1.20 


178 


A586_3 


07 


32:04.5 


+31:39:09.6 


0.41 


1.53 


222 


A586_4 


07 


32:35.4 


+31:35:35.5 


1.03 


-0.47 


227 


A586_5 


07 


32:21.2 


+31:41:26.3 


7.44 


0.46 


232 


A586_6 


07 


32:42.7 


+31:38:37.1 


0.53 


0.16 


293 


A611_0 


08 


01:07.0 


+36:02:18.9 


0.32 


-0.27 


108 


A611_l 


08 


00:52.6 


+36:06:14.2 


0.44 


2.18 


199 


A611_2 


08 


01:17.0 


+36:04:27.8 


0.5 


-0.71 


229 


A621_3 


08 


11:12.8 


+70:02:27.2 


7.18 


1.34 


25 


A621.4 


08 


11:19.3 


+70:00:48.4 


0.6 


-1.31 


127 


A621.5 


08 


11:35.2 


+70:04:25.6 


0.16 


0.13 


166 


A621.6 


08 


10:38.0 


+70:04:09.3 


0.09 


0.01 


181 


A781_0 


09 


20:24.7 


+30:31:49.9 


0.24 


-0.04 


4 


A781_l 


09 


20:23.3 


+30:29:49.3 


8.97 


0.97 


126 


A781.2 


09 


20:08.4 


+30:32:15.8 


1.49 


-0.18 


213 


A781_3 


09 


20:14.0 


+30:28:60.0 


2.12 


0.63 


223 


A990_0 


10 


23:47.3 


+49:11:25.5 


0.44 


-0.21 


208 


A990_l 


10 


24:02.1 


+49:06:51.8 


2.78 


2.18 


239 


A1413_0 


11 


55:15.4 


+23:23:59.4 


0.47 


1.04 


59 


A1413_l 


11 


55:08.8 


+23:26:16.6 


3.1 


0.98 


222 


A1423_2 


11 


57:17.1 


+33:36:30.6 


0.54 


-0.19 


63 


A1423.3 


11 


57:28.5 


+33:35:31.0 


0.26 


2.16 


134 


A1423_4 


11 


57:19.7 


+33:39:58.3 


0.39 


0.73 


171 


A1423.5 


11 


57:35.2 


+33:37:21.8 


0.19 


0.58 


176 


A1423_6 


11 


57:40.5 


+33:35:10.1 


0.18 


-0.40 


270 


A1423_7 


11 


57:20.5 


+33:41:57.8 


0.73 


-0.22 


290 


A1423.8 


11 


57:39.0 


+33:34:03.3 


0.24 


1.47 


290 


A1704_0 


13 


14:02.3 


+64:38:29.5 


0.2 


-0.03 


273 


A1704.1 


13 


14:52.6 


+64:37:59.9 


0.81 


0.86 


287 


A1758a_0 


13 


32:53.3 


+50:31:40.6 


7.08 


0.5 


54 


A1758a_l 


13 


32:38.6 


+50:33:37.7 


0.77 


0.36 


150 


A1758a_2 


13 


33:02.2 


+50:29:26.4 


1.43 


0.28 


190 


A1758a_3 


13 


32:39.6 


+50:34:31.1 


0.3 


1.45 


192 


A1758a_4 


13 


32:41.5 


+50:26:47.7 


0.46 


-0.76 


294 


A1758b_0 


13 


32:33.1 


+50:22:35.1 


0.23 


0.11 


90 


A1758b_l 


13 


32:41.5 


+50:26:47.7 


0.51 


-0.68 


196 


A2009_0 


15 


00:19.7 


+21:22:12.6 


1.85 


3.14 


58 


A2009_l 


15 


00:28.6 


+21:22:45.8 


0.18 


0.66 


133 


A2009.2 


15 


00:19.6 


+21:22:11.3 


1.97 


2.64 


77 


A2009_3 


15 


00:28.6 


+21:22:45.7 


0.18 


0.66 


135 
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Table 5: Source properties for detected sources within 5' of the SZ mean central 
position. The number next to each cluster name denotes the source number; this 
label is used in the plots showing the marginalized posterior distributions for the 
source fluxes. SoMcA is the McADAM-derived best-fit source flux at 16 GHz. a 
is the source spectral index estimated by McAdam and centred at the McADAM- 
derived mean frequency and the last column contains the distance between the 
cluster SZ centroid (as determined by McAdam) and the source. 



Name 


Right ascension 


Declination 


SoMcA 


a 


Distance from SZ centroid 




(hh:mm:ss , J2000) 


(°:':" , J2000) 


mJy 






A2111_0 


15:39:30.1 


+34:29:05.5 


0.5 


0.68 


222 


A2111_l 


15:39:56.7 


+34:29:31.8 


0.81 


-1.47 


297 


A2146_0 


15:56:04.2 


+66:22:13.0 


5.94 


0.55 


43 


A2146_l 


15:56:14.0 


+66:20:53.5 


1.82 


1.03 


59 


A2146_2 


15:56:15.4 


+66:22:44.5 


0.15 


0.34 


89 


A2146_3 


15:55:57.4 


+66:20:03.1 


1.67 


-0.22 


106 


A2146_4 


15:56:27.1 


+66:19:43.8 


0.1 


0.64 


164 


A2146_5 


15:55:25.7 


+66:22:04.0 


0.48 


-0.22 


249 


A2218_0 


16:35:47.4 


+66:14:46.1 


2.86 


0.07 


100 


A2218_l 


16:35:21.8 


+66:13:20.6 


5.99 


0.23 


141 


A2218_2 


16:36:15.6 


+66:14:24.0 


1.77 


0.72 


200 


A2409_0 


22:00:39.7 


+20:58:55.0 


0.75 


1.9 


241 


A2409_l 


22:01:11.2 


+20:54:56.8 


3.12 


0.1 


275 


RXJ0142+2131_0 


01:42:09.2 


+21:33:23.4 


1.09 


0.7 


117 


RXJ0142+2131_1 


01:42:11.0 


+21:29:45.3 


1.16 


1.52 


156 


RXJ0142+2131_2 


01:42:23.3 


+21:30:46.7 


0.3 


0.03 


273 


RXJ1720+2638_0 


17:20:10.0 


+26:37:29.7 


6.92 


1.24 


46 


RXJ1720+2638_1 


17:20:01.2 


+26:36:32.3 


2.05 


0.57 


105 


RXJ1720+2638_2 


17:19:58.4 


+26:34:19.6 


1.22 


1.46 


203 


RXJ1720+2638_3 


17:20:25.5 


+26:37:57.2 


0.88 


0.89 


234 


Zw0857.9+2107_0 


09:00:36.9 


+20:53:41.4 


1.22 


0.31 


102 


Zw0857.9+2107_1 


09:00:55.5 


+20:57:21.2 


0.96 


1.37 


259 


Zw0857.9+2107_2 


09:00:52.8 


+20:58:36.5 


5.57 


0.09 


274 


Zwl454.8+2233_0 


14:57:14.8 


+22:20:34.2 


1.64 


0.28 


14 


Zwl454.8+2233_1 


14:57:08.2 


+22:20:08.6 


1.55 


1.89 


108 


Zwl454.8+2233_2 


14:57:10.6 


+22:18:45.6 


1.49 


0.94 


137 


Zwl454.8+2233_3 


14:56:58.9 


+22:18:49.6 


8.36 


0.17 


258 


Zwl454.8+2233_4 


14:57:04.3 


+22:24:11.9 


0.83 


-0.63 


260 


Zwl454.8+2233_5 


14:57:24.8 


+22:24:52.6 


0.13 


-0.25 


281 


Zwl454.8+2233_6 


14:57:35.7 


+22:19:46.8 


1.04 


1.67 


285 
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Table 6. Mean and 68%-confidence uncertainties for some McADAM-derived large-scale cluster parameters. 



9 



Cluster name 


M T (r 2 oo) 


M T (r 500 ) 


M g (r 2 oo) 


^g(>50o) 


^200 


^500 


Tami 


Ffooo) 


^Osoo) 




xlO u h-^M Q 


xl0 14 h-^M e 


xl0 13 /z7 n 2 n M o 


xlO 13 /*7 o 2 n M 




xlO -1 /^ Mpc 


keV 


xlO 3 arcmur 


xlO 3 arcmnr 


AjoO 


< 1 _i_ i 1 

J.1 ± Z.l 


O 1 _i_ Ci Q 

z.i ± u.y 


A Q _i_ 1 7 
4. J ± 1. / 


z.o ± U. / 


l.Z ± u.z 


0.0 ± l.U 


< 1 _i_ 1 A 
J.Z ±1.4 


3 6 +1 ° 
J.u_ 21 


1 1 _i_ 1 A 
Z. / ± 1 .4 


A611 


IS 


2.0 ±0.5 


3.5 ± 0.6 


2.8 ± 0.3 


1.1 ± 0.1 


6.3 ± 0.5 


4.5 ± 0.6 


2.2 ± 0.5 


2.1 ± 0.4 


A621 


4 8 

^ -1.8 


1.4 ± 0.9 


4.1 ± 1.0 


1.5 ± 0.8 




c q+0.2 

-^-oii 


5.0 ± 1.2 


3 I 

J -1.6 


1.9 ± 1.0 


A773 


3.6 ± 1.2 


1.7 ±0.6 


3 1 +L0 
^ -0.9 


2.1 ±0.4 


1.1 ±0.1 


6.0 ± 0.7 


4 1 +0-9 
-1.0 


1.9 + 0.9 


1 6 +0 - 6 


A/ol 


4.1 ± U.o 


o n _i_ n < 
Z.U ± U.j 


i a _i_ n a 
J.O ± U.o 


Z.y ± U.4 


1 1 _i_ n 1 
1.1 ± U.l 


i _i_ n < 
0.3 ± U.j 


/i < _i_ n 
4. J ± U.O 


^ _i_ (\ a 
Z.j ± u.o 


_i_ 1 r\ < 
Z.Z ± +U.J 


A990 


2 +0 - 4 
z,u -0.1 


1.1 ±0.2 


1.8 ±0.3 


1.6 ±0.2 


0.9 ±0.1 


5.5 ±0.3 


2.8 ±0.3 


0.7 ± 0.2 


0.7 ±0.15 


A1413 


4.0 ± 1.0 


1.9 ±0.6 


3.5 ±0.8 


2.7 ± 0.4 


1.1 ±0.1 


6.6 ± 0.6 


4.4 ± 0.7 


9 9+0.7 
z - -0.8 


2.1 ±0.6 


A1423 


2.2 ±0.8 


1.1 ±0.4 


1.9 ±0.7 


1.5 ±0.4 


0.9 ±0.1 


5.3 ±0.6 


3 +0 - 8 
j. u_ Q7 


0.9 ± 0.5 


0.8 ± 0.4 


A1758a 


4 1+0.7 
^ -0.8 


2.5 ± 0.4 


3.6 ±0.5 


3.4 ± 0.4 


1.1 ±0.1 


6.8 ± 0.4 


4.5 ± 0.5 


2.3 ± 0.5 


2.3 ± 0.4 


A1758b 


4.4 ± 1.9 


2.2 ± 1.0 


o 7+1.6 
^ -1.5 

3.9 ±0.7 


2.2 ± 0.5 


1.1 ±0.2 


6.4 +L1 
6 5 + ° :4> 

° -0.6 


4.6 ± 1.4 


2.7 ± 1.7 


1.6 ±0.6 


A2009 


4.6 + 1.5 


2 +0 - 2 

Z -0.6 


2 4+0.2 
Z -0.3 


1.2 ±0.1 


4 8 +0 - 4 

^ -0.6 


Z -0.7 


2.2 ± 0.9 


A2111 


4.2 ± 0.9 


1.8 ±0.5 


3.6 ±0.7 


2.5 ± 0.3 


1.1 ±0.1 


6.2 ± 0.6 


4.6 ± 0.6 


2.2 ± 0.5 


A2146 


5.0 ±0.7 


2.7 ± 0.5 


4.4 ± 0.5 


3.7 ± 0.4 


1.2 ±0.1 


7.1 ±0.5 


5.2 ±0.5 


3.2 ±0.5 


3.1 ±0.4 


A2218 


6.1 ±0.9 


2.7 ± 0.6 


c 4+O.6 

3 l +d - 9 


4.3 ± 0.4 


1.3 ±0.1 


7.3 ±0.5 


5.9 ±0.6 


4.5 ± 0.7 


4.4 ± 0.7 


RXJ0142+2131 


3 7 +L1 


1.7 ±0.6 


2.1 ±0.4 


1.0 ±0.1 


c q+0.8 


4.2 ± 0.9 


2 +0 - 8 
z.u_ 09 


1 5 +0 - 6 

1 -0.5 


RXJ1720+2638 


2.0 ± 0.4 


1.2 ±0.2 


1.7 ±0.3 


1.6 ±0.3 


0.9 ±0.1 


5.6 ± 0.4 


2.8 ± 0.4 


0.7 ± 0.2 


0.7 ± 0.2 
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6.1 Comparison with gNFW parameterizations 

The adequacy of different profiles, such a s the ft, Navarro Frenk and 
White (NFW), generalized NFW (gNFW jNayarro. Frenk. & W hite 
1996 ) and other hybrid profiles (e.g., [M roczkowskil l201ll 



Allison et al amaie, Hobson, & Grainge 201 1) is still very 



much under debate. We attempt to illustrate the impact that the 
choice of some of these profiles may have on the parameter esti- 
mates by comparing the results obtained from five gNFW param- 
eterizations and from our p parameterization (see Sec. [5} for two 
clusters: Abell 611 (see Sec.[63} and Abell 2111 (see Sec. 157131 . 

For this analysis we sample from the cluster position param- 
eters (x c ,y c ), 6 S = r s /D A , and Y e = Y ioi /D 2 A . r s is the scale radius, 
D A the angular diameter distance and Y tot = Y svh (5r 50 q), where F sph 
is the integrated Compton y parameter within 5r 50 o (lArnaud et al.l 
(12010 ) take 5r 50 o as the radius where the pressure profile flattens). 
Assuming a spherical geometry, 7 sph is calculated by integrating the 
plasma pressure within a spherical volume of radius r: 



^sph(r) = 



—4 

m e c z Jo 



P Q (r')4nr' 2 &v\ 



(6) 



where cr T is the Thomson scattering cross-section, m e is the elec- 
tron mass, c is the speed of light and P e (r) is the electron pressure 
at radius r. The following priors were given for the sampling pa- 
rameters: an exponential prior between 1.3' and 45' for 6 S and a 
power law prior between 0.0005 and 0.2 arcmin 2 for F sph /Z) 2 , with 
a power law index of 1.6 . We note that for the purposes of this ex- 
ercise - to show the different degeneracies for different, plausible 
sets of gNFW-profile parameters - such wide priors are acceptable; 
naturally, where appropriate, the prior ranges can be refined. We 
choose to use a pressure profile for this parameterization since self- 
similarity has been shown to hold best for this q uantity and pressur e 
profiles have low cluster-to-cluster scatter fe.g. jNagai et alJ2007h . 
The gNFW profile is given by 



Table 7. Parameters for the gNFW pressure p rofile. The param eters 
for gNFW^ and gNFWy have been taken from lArnaud et al . 2010 and 
Mroczkowski et al. 2009, respectively. The values in Mroczkowsk^UilJ 
2009 are the corrected values for the results published by Nag ai et al.l2007l . 



Profile Label 


a 


b 


c 


C500 


gNFWi 


1.37 


5.49 


0.035 


2.16 


gNFW 2 


0.33 


5.49 


0.065 


0.17 


gNFW 3 


2.01 


5.49 


0.860 


1.37 


gNFW A 


1.0620 


5.4807 


0.3292 


1.156 


gNFWjv 


0.9 


5.0 


0.4 


1.3 




Figure 2. 2-D marginalized distribution for y S ph( r 50o) against r^oo obtained 
using the /3-based cluster parameterization without any data. 
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P ei , the overall normalisation coefficient of the pressure profile, is 
calculated by computing Eq. © for 5r 50 o; once P e i has been found, 
y can be obtained. The shape of the gNFW profile is governed by c 
in the inner cluster regions (r « r s ), by a at intermediate radii (r « 
r s ), and by b on the cluster outskirts (r > r s ). These parameters, 
together with the concentration param eter, c 50 n, are fixed i n most 
analyses to some best-fit values (e.g.. iMroczkowskil 1201 lh . With 

r s can be expressed in terms of r 5 oo : r s = rsop/csoo, w hich 
is a commo n reparameterization (see e.g.. lArnaud et al.lkoiol and 
iNagai et alJ200% . 

We ran our analysis using a gNFW profile with parameters 
defined by Nagai et al. as gNFW N , another defined by Arnaud et al. 
as the 'universal' profile gNFW A , and three other combinations for 
the slope parameters and c 50 o that were found to provide the best fit 
for some clusters in Arnaud et al.; gNFW l9 gNFW 2 and gNFwJl 
The gNFW parameters for our five choices are given in Tab. [7] 

The 2D-marginalized posterior distributions of Fsphfooo) 
against r 50 o obtained for each of the five parameterizations, as well 
as the p parameterization from Sec. [5] for Abell 611 and Abell 2111 
are shown in Fig.Q]. We test for possible biases in our results from 
the choice of priors by running the analysis without data; the results 



8 The parameters for the three other gNFW profiles all lie within 3<x of 
the average value of each parameter obtained using the cluster sample in 
Arnaud et al. . 



indicate the constraints imposed by our priors. We find no evidence 
for significant biases, as shown in Fig. [2] 

When the shape parameters of the f3 profile are fitted to the SZ 
data instead of being set to the X-ray value (typically derived for 
data sensitive to smaller scales than AMI data) we find the mean 
values for F sph (r 5 oo) and r 50 o derived from the p analysis to be con- 
sistent (within 1-2<t) with those from gNFW A and gNFW N - the 
averaged gNFW profiles. For these two clusters we find all gNFW 
parameterizations yield lower values for F sph (r 5 oo) than for the p 
analysis; this is not the case for r 50 o, for which no systematic dif- 
ference is seen. The constraints on F sph (r 5 oo) are similar for most of 
the gNFW models (with the exception of gNFW 2 ) and the p model, 
while those for r 50 o appear to be tighter for the p model. One strik- 
ing difference between the two types of parameterizations is the 
shape and orientation of the F sph - r 50 o degeneracy. 

The resolution and limited spatial dynamic range of the AMI 
data do not allow profile selection to be made robustly, as indicated 
by the small difference in evidence values between the different pa- 
rameterizations (Tab. [8j. Hence, our p parameterization provides 
a comparable fit to that of the commonly used, averaged gNFW 
profiles, gNFW A , and gNFW N . It is clear from Fig.[T|that the dis- 
tribution for the r sph (r 5 oo) - r 50 o degeneracies is very sensitive to 
the choice for the slope parameters (and c 50 o for gNFW). Cluster 
parameters for a cluster with a profile described by e.g., a gNFW 2 
recovered using a gNFW A paramaterization will be biased. 
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Figure 1. 2-D marginalized distributions for y S ph(^50o) against r^oo obtained using the /3-based cluster parameterization and five gNFW-based cluster param- 
eterizations with slope parameters and C500 given in Tab. [7] The crosses denote the McADAM-derived mean values. The results are for Abell 611 (left) and 
Abell 2111 (right). The blue filled ellipses show the results of the [3 parameterization. 



Table 8. Log e evidences for five cluster parameterizations applied to 
Abell 2111 and Abell 611 



Abell 2111 Abell 611 



Profile Label 






gNFWi 


23198.88 


21114.08 


gNFW 2 


23199.51 


21114.40 


gNFW 3 


23198.64 


21114.74 


gNFW A 


23198.94 


21114.50 


gNFW^v 


23198.76 


21114.65 




23194.92 


21112.05 



6.2 Abell 586 

Results for Abell 586 are given in Figs. |3]and|4] This cluster has 
a complex source environment, with 7 sources within 5' from the 
cluster SZ centroid, which include two radio sources of « 260 
and 744 yt/Jy at 0.5' and 4' from the pointing centre. After source- 
subtraction there are only « lcr residuals left on the map. Uncer- 
tainties in the source fluxes are carried through into the posterior 
distributions for the cluster parameters. From Fig. @] it can be seen 
that there is no strong degeneracy between the source flux densities 
and the cluster mass. 

Abell 5 86 h as been studie d extensively in the X-ray band (e.g., 
lAllenl 120001 and IWhite] |2000). A recent analysis of the tempera- 
ture profile dCvprianoll2005l) shows how the temperature falls from 
« 9 keV at the cluster centre to « 5.5 keV at a radius « 280" . Cypri- 
ano et al. have used the Gemini Multi- Object Spectrograph together 
with X-ray data taken from the Chandra archive to measure the 
properties of Abell 586. They compare mass estimates derived from 
the velocity distribution and from the X-ray temperature profile and 
find that both give very similar results, M g « 0.48 x 10 14 M o (for 
hjo = 1) within Mpc. They suggest that the cluster is spheri- 

cal and relaxed with no recent merg ers. It is less clear w hether this 
cluster has a cool core or not, with lAUen & Fabian ( 1998) reporting 



its existence and Marrone et al ] (EoT 1 ) saying otherwise. The peak 
X-ray and SZ emissions are consistent with each other and the AMI 
SZ decrement shows some signs of being extended towards the S W 
(Fig.[3]B and C); there are no contaminating sources in the vicinity 
of this SZ-'tail'. 

The SZ effect fr om Abell 586 h as previously been observed 
with OVRO/BIMA bv lLaRoaue et al] <2006h and lBonamente et alJ 
(2006). LaRoque et al. apply an isothermal /?-model to SZ and 
Chandra X-ray observations and find M g (r 25 oo) = 2.49 + 0.32 X 
10 13 M o and M g (r 25 oo) = 2.26^ x 10 13 M o , respectively (us- 
ing h 10 = 1 and excising the inner lOOkpc from the X-ray data). 
In addition, they determine an X-ray spectroscopic temperature of 
the cluster gas of « 6.35 keV be t ween a radius of lOOkpc and 
r 2 5oo- In comparison, lOkabe et al.l (120 101) use Subaru to calculate 
the cluster mass from weak lensing by applying a Navarro, Frenk 
& White (NFW: lNavarro. Frenk. & Whitell 19961) profile. They find 
M T (r 25 oo) = 2.41^ x 10 14 M o and M T (r 500 ) = 4.74!}$ x 1O 14 M 
(using h 12 = 1). In this work, we find M T (r 500 ) = 3.0+1.3xl0 14 M o , 
where r 500 = 0.94+0. 14 Mpc and h 10 = 1 . Note that the fluxes of the 
radio sources S 1 and S 5 are degenerate in our analysis of Abell 586 
(see Fig. [4}; this is because these sources are separated by only 66" 
and their individual fluxes cannot be disentangled in the analysis of 
the AMI SA data. 



6.3 Abell 611 

Results for Abell 621 are presented in Fig. [5] Our methodol- 
ogy is able to model the radio sources + cluster enviroment 
well, as demonstrated by the good constraints on the mass and 
other parameters and the lack of degeneracies between the sources 
closest to the cluster and the cluster mass (Fig. [5] D, E and 
F). We do not expect any significant contamination from radio 
sources nor f r om ex tended emission since GMRT observations by 
IVenturi et all (2008) found no evidence for a radio halo associ- 
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Figure 3. Results for Abell 586. Panels A and B show the SA map before and after source subtraction, respectively; a 0.6 kA taper has been applied in B. The 
box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] Panel C shows the smoothed Chandra X-ray map overlaid with 
contours from B. D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. The v-axis for the 1-d 
marginals is the probability density and for all the posterior distributions plots in this paper h refers to h\oo. In D Mj is given in units of h^ x 10 14 M o and 
f g in /z~q ; both parameters are estimated within r2oo- in E M g is in units of h~Q Q M Q , r in h~^ Mpc and T in KeV. The slight rise in the distribution for r2oo at 
large r is a result of a binning artifact and, in fact, this distribution does tail off smoothly, as expected. 

© 0000 RAS, MNRAS 000, 000-000 



SZ study of 19 LoCuSS galaxy clusters with AMI out to the virial radius 13 




5 10 15 7 7.5 
M T ( r 20 ) S 5 



Figure 4. 1 and 2-D marginalized posterior distributions for the flux den- 
sities, in Jys, of sources detected within 5' of the SZ centroid of Abell 586 
(see Tab.[5} and Afxfaoo), in units of h^ Q M Q x 10 14 . 



ated with Abell 611 at 610 MHz. The decrement on the source- 
subtracted maps appears to be circular, in agreement with the 
X-ray surface brightness from the Chandra achive data shown 
in Fig. [5] C, which also appears to be smooth and whose peak 
is close to the position of the brightest cluster galaxy and the 
SZ peak. These facts might be taken to imply th e cluster is re- 
laxed but, it does not seem to have a cool core jMarrone et al.l 
2011). Abe ll 611 has also previo usly been observed in the SZ at 
15 GHz bv [Grainger et all (120021). AMI Consortium: Zw art et al. 
(2010) and lAMI Consortium: Hurlev-Walker et all d201lh. an d at 
30 GHz by [ Bonament e et all <2004h . lBonamente et alf (2006) and 
iLaRoaue et al.U2006h . 

From the analysis in iDonnarumma et al ] (Soul the clus- 
ter mass was estimated to be 9.32-1 l.llxl0 14 M o (within a ra- 
dius of 1.8+0.5 Mpc) by fitting different cluster models to X- 
ray data and between 4.01-6.32xl0 14 M o (within a radius of 
1.5+0.2 Mpc) when fitting different models to the lensing data; 
all estimates use h 10 = 1. Several other analys es of Chandra 
data produce comparable mass estimates (e.g.. ISchmidt & Allen 
l200llMorandi. Ettori & Mosca rdini 2007, Morandi & Ettori1 l2007l 



and lSanderson et al.ll2009l) . lRomano et al. ( 2010h perform a weak- 
lensing analysis of Abell 611 using data from the Large Binocu- 
lar Telescope; with an NFW profile they estimate M T (r 2 oo) = 4- 
7xl0 14 M o for/* 70 = 1. 

These are in agreement with the values obtained from Sub- 
aru weak lensing observations by Okabe et al.. AMI Consortium: 
Hurley-Walker et al. estimate the total mass for this system within 
r 2 oo- Using lensing data they find it is 4.7 + 1.2 x 10 14 /z^M o and 
using AMI SZ data they find it is 6.0 + 1.9 X 10 14 ^ 1 M o . We find 
M T (r 20 o) = 5.7 + 1.1 x 10 14 M o , where r 20 o = 1.6 + 0.1 and h 10 = 1; 
this value is significantly smaller than the result given in AMI Con- 
sortium: Zwart et al. (2010); this is due to their mass measure- 
ments being biased high, as they said, and is further discussed in 
lAMI Consortium: Olamaie et all (1201 ll) . 



6.4 Abell 621 

Fig-Hlcontains our results for Abell 621. Out of the 13 radio sources 
detected on the LA raster for Abell 621, three lie near the edge of 
the cluster decrement in the source- subtracted map and one, which 
has a flux density « 7 mJy, is coincident with the best-fit cluster 
position, as indicated by the box in Fig. [6] A. However, whatever 
reasonable source subtraction we try makes almost no difference to 
the inferred cluster mass. The ROSAT HRI X-ray image presented 
in Fig.[6]C appears to be quite uniform and circular and the offset 
between the X-ray and SZ cluster centroids is small. We find the 
cluster mass to be M T (r 2 oo) = 4.8^} gXlO 14 /i~ o 1 o M from our analy- 
sis; at 6cr, this is one of our less significant detections. 

The data for the probability distributions in Fig.[6]E have been 
binned relatively finely to avoid misleading features, in particu- 
lar towards the lower limits of our plots. As a result, the noise 
in these bins is higher, which makes the distributions appear less 
smooth. For some combinations of cluster parameters, there is 
nowhere in the cluster density estimation that the density of the 
gas reaches 500p crit . In these cases, where there is no physical so- 
lution for r 5 oo, we set r 50 o = 0. This leads to sharp, meaningless 
peaks at small radius in the distributions for some cluster parame- 
ters at rspp (Fig.[6lE). These feature s have also been discussed in 
AMI Consortium: Zwart et all j2010h . 



6.5 Abell 773 

Results for Abell 773 are shown in Fig. [7] Abell 773 has few 
associated radio sources, all of which are > 10' away from the 
pointing centre, weak (< 3 mJy), and are subtracted well from our 
data (Fig|7]B). We do not find any ev idence for extended positive 
emiss ion in our maps. Observations bv lGiovannini. Tordi. & FerettH 
(119991) revealed the presence of a radio halo with a luminosity of 
2.8 X 10 24 WHz" 1 at 1.4 GHz; th is result has been confirmed with 
the VLA bv iGovoni et al .1 d200 lh . Given the typical steep spectral 
index of radio halos, we do not expect our SZ signal to be affected 
at 16 GHz. 

Our observations clearly show the SZ image is extended along 
the NW-SE direction, contrary to the X-ray image from Chandra 
observations, which appears to be elongated in an approximately 
perpendicular direction. As might be expected fr om a disturbed 
system, Abell 773 appears to not have a cool core dAUen & F abian 
Q991L 

Barrena et al ] (|2007b present a comprehensive study of 
Abell 773 from the Telescopio Nazionale Galileo (TNG) telescope 
and X-ray data from the Chandra data archive. They find two 
peaks in the velocity distribution of the cluster members which are 
separated by 2' along the E-W direction. Two peaks can also be 
seen in the X-ray, although these are along the NE-SW direction. 
Barrena et al. estimate the virial mass of the main cluster to be 
M T (r vir ) = 1.0 - 2.5 x \0 15 h^M Q and M T = 1.2-2JxlO l5 h^M Q 
for the entire system, using the virial theorem, dispersion veloc- 
ity measurements and a galaxy King-like distribution. Assuming 
an NFW profile they estimate the mass for the system to be M T (< 
r = l/^Mpc) = 5.9 - 11.1 x 10 14 h ^Mft. A further analysis of 
Chandra data bv lGovoni et al.l J2004I) yielded a mean temperature 
of 7.5+0.8 keV within a radius of 800 kpc (h 10 = 1). Another X-ray 
study of this cluster by Zhan g et aD d2008h using XMM-Newton 
found M T (r 5 oo) = 8.3 +2.5 xl0 14 M o assuming isothermality, spher- 
ical symmetry and h 10 = 1. 

The SZ effect associated with Ab e ll 773 has been ob- 
served several times ( Grain ge et all Il993l ICarlstrom et all ll 9921 . 
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Figure 5. Results for Abell 61 1. Panels A and B show the SA map before and after source subtraction, respectively; a 0.6 kA taper has been applied to B. The 
box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] C shows the smoothed Chandra X-ray map overlaid with contours 
from B. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F shows the 1 and 2-D 
marginalized posterior distributions for source flux densities, in Jys, (within 5' of the cluster SZ centroid, see Tab. [5} and Mjfooo), in h^L x 10 14 M o . In panel 
D Mj is given in units of h~L x 10 14 M o and f g in h~L; both parameters are estimated within r2oo- In E M g is in units of /z"iM , r in h~LMpc and T in KeV. 
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Figure 6. Results for Abell 621. Panels A and B show the SA map before and after source-subtraction, respectively; a 0.6 kA taper has been applied to B. The 
box in panels A and B indicates the cluster SZ centroid, the other symbols are in Tab.@] The smoothed ROSATHRI X-ray map overlaid with contours from B 
is given in panel C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F shows the 1 
and 2-D marginalized posterior distributions for source flux densities (in Jys) within 5' of the cluster SZ centroid (see Tab.|5} and Afxfooo) in h^ Q x 10 14 M o . 
In D Mj is given in units of h^ x 10 14 M© and f g in h^ Q ; both parameters are estimated within r2oo- In E M g is in units of h^ Q M Q , r in /z^Mpc and T in 
KeV. 

© 0000 RAS, MNRAS 000, 000-000 



16 AMI Consortium: Rodriguez-Gonzdlvez et al. 



ISaunders et al.l2003LlBonamente et alkood LaRoque et al. 2006). 
Most recently, AMI Consortium: Zwart et al. (2010) observed the 
cluster and found a cluster mass of M T (r 2 oo) = 1.9^Q4Xl0 15 M o us- 
ing &ioo = 1; however, their M T estimates are biased high, as they 
say, and we find M T (r 2 oo) = 3.6 ± 1.2 x 10 14 /*i oM . 



6.6 Abell 781 

Fig. [8] contains our results for Abell 781. It is evident from inspec- 
tion of Figs. [8] A and F and Tab. \5\ that there is strong emission 
from radio sources lying on the decrement. One of the sources with 
a flux density of 9 mJy lies on top of the McAdam best-fit cluster 
position. The difficulty of accurately disentangling the signal con- 
tributions from this source and the cluster is translated into a degen- 
eracy between the source's flux density and the cluster mass: Fig. 
[8]F. No extended emission was detected on the LA maps and, after 
source subtraction, the residual s on the maps are < 2<x (Fig.[8]B). 
Rud nick & Lemmermanl d2009h have found evidence in WENSS 
data at 327 MHz of diffuse emission from a radio galaxy and some 
other unknown source with a flu x within a radius of 500 kpc of 
40 mJy, while Ventu ri et al 1 d2QQ8h estimate diffuse emission at the 
centre to be « 15 - 20mJy using 325-MHz GMRT data. Assum- 
ing a typical st eep spectral in dex for radio halos, in the range of 
1.2-1.4 (e.g., Hanisch 1980), even as far as 16 GHz, we would 
expect to find an « 170/Jy signal around the cluster and < 85 /dy 
at the centre. The GMRT contour map in Venturi et al. identifies 
the relic at a similar location to that of some unsubtracted positive 
emission in our maps at * RA 09:30:00, Dec 30:28:00. 

X-ray observ ations with Chandra and XMM-Newton 
(Seh galet ai]l2008h imply that Abell 781 is a complex cluster 
merger: the main cluster is surrounded by three smaller clusters, 
two to the East of the main cluster and one to the West. Sehgal 
et al. estimate the mass of Abell 781 within r 50 o assuming a NFW 
matter density profile to be 5.2^xl0 14 M o from X-ray data and 
2.7+ 1 y°xl0 14 M o from the Kitt Peak Mayall 4-m telescope lensing 
observations, using hj\ = 1. Further results from XMM-Newton by 
Zhang et al. yield M T (r 50 o) = 4.5+1.3 xl0 14 M o assuming isother- 
mality, spherical symmetry and /z 70 = 1. We obtain M T (r 50 o) = 
2.9+O.6xlO 14 M and M T (r 200 ) = 5.9+1. lxl0 14 M o for h 10 = 1. 



6.7 Abell 990 

Results for Abell 990 are given in Fig. [9] We detected 20 sources 
towards Abell 990. Those detected above 4cr LA within 10' from the 
pointing centre were found to have flux densities < 2.8 mJy, not to 
be extended with respect to the LA synthesized beam (Tab. [5]), and 
none to lie on the SZ decrement, as seen in the source- subtracted 
map (Fig.|9]B). The subtraction has worked well and there are only 
low-level (« 1 — 2cr) residuals. iRudnick & Lem merman (20091) do 
not detect any significant amount of diffuse emission within a ra- 
dius of 500 kpc in 327 MHz WENSS data; given the steep falling 
spectrum associated with this emission, we do not expect it to con- 
taminate our SZ signal. The imaged decrement is fairly circular 
but extended along the NE-SW direction coincident with the distri- 
bution of the X-ray signal. Our spherical cluster model provides a 
good fit and the parameter distributions are tightly constrained. The 
low resolution X-ray map shown in Fig.|9]C provides tentative evi- 
dence that the X-ray emitting cluster gas has a clumpy distribution. 



6.8 Abell 1413 

In Fig. [10] we present results for Abell 1413. It can be seen from 
Figs. [10] A and B that there are two of sources on the decrement 
with flux densities of 0.47 and 3.1 mJy (in Tab. \5}. The brightest 
source in our LA maps has a flux density of 14mJy but, since it 
is 700" from the cluster X-ray centre, it does not contaminate our 
SZ signal. Some residual flux is seen on the source- subtracted SA 
maps; the strongest residuals are not associated with sources in the 
LA data, suggesting they could be exten ded emission resolved out 
from the LA maps iGovoni et al] d2QQ9h find tentative (« 3(f) evi- 
dence in FIRST data at 1.4 GHz for a weak mini halo with a lumi- 
nosity of 1.0 X 10 23 W Hz -1 . The peak signal from this mini halo is 
offset to the East with respect to the central cD galaxy, similarly to 
our SZ peak, which is slightly offset to the SE of the X-ray centroid. 
Abell 1413 does seem to be a relaxed cluster; this is supported by 
the smooth X-ray distribution, the good agreement between the X- 
ray and SZ centroids, the circular appe arance of the projecte d SZ 
signal and the presence of a cool core dAllen & Fa bian 1998b. We 
therefore expect our model to provide a good fit to the AMI data 
towards this cluster. 

Abell 1413 has been observed in the X-r ay by XMM-Newton 
(e.g .. IPratt & Arnaud|l2005h . Chandra (e.g., IVikhlinin et al]l2003 
and|gonamente et al. 2006) and most recently by the Suzaku satel- 
lite dHoshino et al. 2010); SZ images have been made with the Ryle 
Telescope at 15 GHz Jdrainge et al.lll996h and with O VRO/BIMA 
at 30 GHz (LaRoque et al. and lBonamente et alj|20060 . These anal- 
yses indicate that Abell 1413 seems indeed to be a relaxed clus- 
ter with no evidence of recent merging. Different temperature and 
density profiles obtained from X-ray data are in good agreement 
out to half the virial radius. Hoshino et al. measure the variation 
of temperature with radius, finding a temperature of 7.5 keV near 
the centre and of 3.5 keV at r 2 oo; they assume spherical symmetry, 
an NFW density profile and hydrostatic equilibrium to calculate 
M T (r 20 o) = 6.6+2.3xl0 14 /*7 1 M o . Zhang et al. use XMM-Newton 
and find M T (r 500 ) = 5.4 +1.6 xl0 14 M o ; they assume isothermal- 
ity, spherical symmetry and h 7 o = 1. We determine M T (r 2 oo) to be 
5.7 + 1.4 x 10 14 M o for/* 70 = 1. 



6.9 Abell 1423 

Results for Abell 1423 are shown in Fig. Qj] The source en- 
vironment for Abell 1423 is challenging (see Fig. [TT] A) - 23 
sources have been detected within 10' of the X-ray cluster cen- 
troid, of which 4 lie on the decrement, as seen from the source - 
subtracted map. We find no evidence for extended emission, in 
agreement with the la ck of diffuse emission towards t his cluster at 
327 M Hz reported by IRudnick & Lemmermanl d2009h and the re- 
sults in lRossetti et al.1 fcoill) . The sources closest to the cluster all 
have flux densities < 1.3 mJy (Tab. [3 and only small positive resid- 
uals remain after source subtraction. As shown in Fig. [TTlR the flux 
densities for some of the sources close to the cluster centroid man- 
ifest degeneracies with the cluster mass. 

The details on the dynamics of Abell 1423 are largely un- 
known. The lack of strong radio halo emission i s indicative o f a 
system without very significant dynamical activity dBuot e 2001), as 
is the good agreement between the X-ray and SZ emission peak po- 
sitions. On the other hand, the X-ray data in Fig.Qj]C shows signs 
of substru cture and our SZ image is be elongated along the SE-NW 
direction. IS anderson et al.l (12009b find that the logarithmic gradient 
for the gas density profile of Abell 142 3 at 0.04r5oo is a » -0 9 8 - a 
key signature of cooling core clusters (Vikhlinin e7aLl(l2006h sug- 
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Figure 7. Results for Abell 773. Panels A and B show the SA map before and after source subtraction, respectively; a 0.6 kA taper has been applied to B. The 
box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] The smoothed Chandra X-ray map overlaid with contours from B 
is given in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. In panel D 
Mt is given in units of h^ Q x 10 14 M o and f g in /z~J ; both parameters are estimated within r2oo- In E M g is in units of h^ M Q , r in /z^Mpc and T in KeV. 
No plots of the degeneracy between cluster mass and source flux densities are shown since all detected sources are > 5' from the cluster SZ centroid and thus 
should not have a strong impact on the marginalized distribution for the cluster mass. 
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Figure 8. Results for Abell 781. Panels A and B show the SA map before and after source-subtraction, respectively; a 0.6 kA taper has been applied to B. The 
box in panels A and B indicates the cluster SZ centroid, other symbols are in Tab. [4] The smoothed Chandra X-ray map overlaid with contours from B is given 
in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F shows the 1 and 2-D 
marginalized posterior distributions for source flux densities (in Jys) within 5' of the cluster SZ centroid (see Tab. [5) and Mx(r2oo) (in h^ Q x 10 14 M o ). In D 
Mt is given in units of h^ x 10 14 M© and f g in h^ ; both parameters are estimated within r2oo- In E M g is in units of h^ M o , r in /z^Mpc and T in KeV. 
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Figure 9. Results for Abell 990. Panels A and B show the SA map before and after source-subtraction, respectively; a 0.6 kA taper has been applied to B. The 
box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab.|4] The smoothed ROSAT HRI X-ray map overlaid with contours from 
B is shown in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F shows 
the 1 and 2-D marginalized posterior distributions for source flux densities in Jy (see Tab. [5} and Mt(>2oo) (in xl0 14 M o ). In panel D Mj is given in units of 
h^Q Q x 10 14 M© and f g in h^ ; both parameters are estimated within r2oo- In E M g is in units of h^ M Q , r in /z^Mpc and T in KeV. 
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Abell 1413 




Figure 10. Results for Abell 1413. Panels A and B show the SA map before and after source-subtraction, respectively; a 0.6 kA taper has been applied to 
B. The box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] The smoothed Chandra X-ray map overlaid with contours 
from B in presented in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. 
F shows the 1 and 2-D marginalized posterior distributions for source flux densities (in Jys) within 5' of the cluster SZ centroid (see Tab. [5) and Mx(r2oo) (in 
h~Q Q x 10 14 M o ). In panel D Mj is given in units of h~Q Q x 10 14 M o and f g in h~^ Q \ both parameters are estimated within T2qq. In E M g is in units of H'^Mq, r 
in /z^Mpc and T in KeV. 
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gest a < -0.7 for strong cooling flows). In their study clusters with 
small offsets at r 50 o between the X-ray and the Brightest Cluster 
Galaxy (BCG) are tightly correlated with large, negative spectral 
indices, an indication that the strength of cooling cores tends to 
drop in more disturbed systems, but Abell 1423 is an unsual outlier 
in this trend with a small offset and a steep a. 



6.10 Abell 1704 

Abell 1704 has been observed with ROSAT HRI and PSPC 
teizza el al J 1993) . These observations show a shift in position be- 
tween the peak emission and the cluster centroid and distinct signs 
of elongations in the gas distribution. Further analy sis of X-ray ob- 
serva ti ons suggest the presen ce of a cooling flow dAUen & Fabianl 
Il998h . lCarlstrom ^taiT(ll992h attempted to detect an SZ effect us- 
ing the OVRO array at 30 GHz towards this cluster but found no 
convincing SZ signal. 

The NVSS map at 1.4 GHz shows complex, extended emis- 
sion (Fig. [12]). These features are detected on our SA maps but a 
significant portion of the emission is resolved out on our LA maps 
(see Tab. [5] for more details on these sources). Our model is not 
sophisticated enough to deal properly with extended structure and 
significant residual emission can be seen in the source- subtracted 
SA map, Fig. [12] Consequently, we are not able to convincingly 
detect an SZ effect towards Abell 1704. 



6.11 Abell 1758 

Results for Abelll758a and b are given in Fig. [l3}Fig [15] It 
is clear from Fig. [13] A, B and C that Abell 1758 is a com- 
plex system comprising two gravitationally-bou nd main clusters, 
Abell 1758a and Abell 1758 b, separated by 8' dRizza el al.lfl998l 
and iDavid & Kempnej|2004h . David & Kempner find no conclu- 
sive evidence for interaction between these two main clusters, yet 
each of them is undergoing major mergers - Abell 1758a between 
two 7-keV clusters and Abell 1758b between two 5-keV clusters; 
since both sets of mergers are between clusters of approximately 
equal mass, provided each of the primary clusters was virial ized 
pre-merger, we might expect the average temperature to be higher 
by some 25% when all the gas mass of the subcluster has merged 
with that of the primary cluster. 

To map the full extent of this system we took raster observa- 
tions with the SA, which are presented in Fig.[l3]B and C. From[T3l 
Cii it can be seen that the SZ signal follows the X-ray emission but 
there seems to be a hint of an SZ signal connecting these two clus- 
ters; note that the clusters have identical redshifts. No connecting 
X-ray signal would be expected and ind eed none is seen. A recent 
analysis of Spitzer/MIPS 24^m data bv lHaines et all d2009h classi- 
fies Abell 1758 as the most active system they have observed at that 
wavelength. They also identify numerous smaller mass peaks and 
filamentary structures, which are likely to indicate the presence of 
infalling galaxy groups, in support of the David & Kempner obser- 
vations. 

For Abell 1758a we obtain M T (r 500 ) = 2.5 +0.4 x \0 1A h m M Q 
and M T (r 20 o) = 4.lg;| x lO u h m M Q . Zhang et al. studied 
Abell 1758a using XMM-Newton and found M T (r 50 o) = 1.1 ± 
0.3 x 10 15 M o ; they assumed isothermality, spherical symmetry and 
h 10 = 1. 




Figure 15. 1 and 2-D marginalized posterior distributions for Mjfooo) and 
sources within 5' from the cluster X-ray centroid for Abell 1758a. Source 
flux densities are given in units of Jys and Mx(r2oo) in units of h~^ M Q x 
10 14 . 

6.12 Abell 2009 

Results for Abell 2009 are given in Fig. [16] Eighteen sources were 
detected above 4cr LA in our LA maps. Given that all of the sources, 
except one, are further away than one arcminute from the point- 
ing centre and have flux densities < 2mJy, the source environ- 
ment should not significantly contaminate the SZ signal on the SA 
maps. The source- subtraction has worked well and there are only 
2cr residuals (Fig.[l6] B); the most prominent residual is likely to be 
associated with some extended emission seen in the SA map before 
source subtraction. 

We fi nd the SZ image is extended in an approximately NS 
direction. Okabe et a D boioh fit an NFW profile to weak lensing 
data from the Subaru/Suprime-Cam and find M T (rn ) - 3.86+J- 93 x 
10 14 /^M o (with h 72 = 1.0). We find M T (r 20 o) = 4.6 + 1.5 x 
lO^Mo. 

The misleading sharp peaks at small radius in the distributions 
for cluster parameters at r 500 (Fig.[l6]F) are discussed in Sec. 16.41 

6.13 Abell 2111 

Results for Abell 2111 are presented in Fig. [T7] The source en- 
vironment in the vicinity of Abell 2111 does not present a prob- 
lem in our analysis: all the sources are located on the edge of the 
decrement or beyond and have flux densities < 3 mJy (Fig. [T7] 
A and Tab, [5}. Some residual flux with a peak surface brightness 
« 700 /Jy beam -1 remains in our source- subtracted map but is suf- 
ficiently far (« 45") that it has a negligible effect on our SZ detec- 
tion (Fig. [T7|B). 

X-ray studies of ROSAT PSPC and HRI data by IWang et al.l 

dl997h reveal Abell 2111 has substructure on small scales but ap- 
pears to be reasonably relaxed on larger scales away from the core. 
Wang et al. identify a main X-ray emitting component and a hot- 
ter subcomponent and conclude that Abell 2111 is most likely to 
be a head-on merger be tween two subclusters; this is supported by 
Henri ksen etail (1999) using ASCA data. A disturbed nature of 
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Figure 11. Results for Abell 1423. Panels A and B show the SA map before and after surce-subtraction, respectively; a 0.6 kA taper has been applied to B. 
The box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] The smoothed Chandra X-ray map overlaid with contours 
from B is presented in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. 
F shows the 1 and 2-D marginalized posterior distributions for source flux densities (in Jys) within 5' of the cluster SZ centroid (see Tab. [5} and Afrfooo) (in 
h^Q Q x 10 14 M o ). In panel D Mj is given in units of h^ Q x 1O 14 M and f g in h^ Q ; both parameters are estimated within r2oo- In E M g is in units of h^ Q M Q , r 



in /z 1( j Mpc and T in KeV. 
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Figure 12. A: source-subtracted SA map produced using a 0.6-k/l taper. The contours increase linearly in units of osa- B: 1.4-GHz NVSS map towards Abell 
1704. 



Abell 2111 might also be indicated by the apparently clumpy X- 
ray emission and X-ray-SZ offset seen in Fig.fTTlC. 

Recent investigations by iRines et aTl J2010h find the virial 
mass for Abell 2111 to be M T (r 100 ) = 4.01 + 0.41 x 10 14 M o us- 
ing hjo = 1.0 from an av erage of 90 member redshifts within 
rioo. iMaughan et al ] d2008h fit a modified version of the stan- 
dard ID isothermal /?-model to Chandra data with h 10 = 1.0 
to compute M g (r 50 o) = 7.44+° o° x 10 13 M o . We obtain a value 
of Af g (r 5 oo) = 2.5 + 0.3 x \O u h^M Q . Previously. ILaRoaue et all 
(2006) fitted an isothermal /5-model to Chandra data (excising the 
r < lOOkpc from the core) and OVRO/BIMA data and found a 
gas mass M g (r 25 oo) = 2.15 + 0.42 x 10 13 M o (for h 70 = 1.0); they 
also found an X-ray spectroscopic temperature of « 8.2 keV. On 
larger scales, at r 20 o, we obtain a lower temperature, 4.6 + 0.6 keV, 
which suggests the average cluster temperature falls with radius. 
Moreover, Henriksen et al. report a radially decreasing temperature 
structure for Abell 2111 and parameterize it by a poly tropic index 



y « 1 .45. On larger scales [AMI Consortium: Hurley-Walker et al.l 
< 201lh estimate M T (r 20 o) = 6.9+ 1.1 X 10 u h~^M Q from lensing data 
and M T (r 20 o) = 6.3 + 2.1 x 10 14 ^ 1 M o from AMI SZ data; they also 
find that a circular geometry is a slightly better fit to the data than an 
elliptical geometry. Our results, M T (r 20 o) = 4.2 + 0.9 x 10 u h~Q Q M Q 
are in very good agreement. 



6.14 Abell 2146 

We have re-analysed the AMI data used in 
AMI Consortium: Rodriguez-Gonzalvez et al] d201Qh with the 
cluster parameterization described in Sec. 15. lj which is slightly 
different to theirs; our results are presented in Fig. 16.141 They 
obtain M g (r 200 ) = 4.9 + 0.5 x \0 l3 h^ M Q and T = 4.5 + 0.5 keV 
while our results give M g (r 20 o) = 4.4 + 0.6 X 10 13 /*~q M o and 
T = 5.2 + 0.5 keV. Given the similarities between the two analyses 
and the fact the same data were used for both, we would indeed 
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Figure 13. Panels A. and B. show the SA map before and after source subtraction (the latter map has had 0.6-kA taper applied to it). The boxes in panels A 
and B indicates SZ centroid for each cluster, for the other symbols see Tab. [4] The maps shown here are primary beam corrected signal-to-noise maps cut off 
at 0.3 of the primary beam. The noise level is « 115//Jy towards the upper cluster (Abell 1758a) and « 130//Jy towards the lower cluster (Abell 1758b). The 
source-subtracted SA maps from B are overlaid with the the Chandra map in Ci and with ROSAT PSPC X-ray map in Cii. D and E show the marginalized 
posterior distributions for sampling and derived parameters, respectively. In panel D, Mj is given in units of h^ Q x 10 14 M o and f g in h~^ Q M Q ; both parameters 
are estimated within r2oo- I n E M g is in units of ^ O M , r in /z^Mpc and T in KeV. 
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Figure 14. Abell 1758b. Left panel: 1 and 2-D marginalized posterior distributions for the cluster sampling parameters. Mj is given in units of h^ x 10 14 M© 
and f g in h^ ; both parameters are estimated within r2oo- Right panel: 1-D marginalized posterior distributions for the cluster derived parameters. M g is in 
units of h~Q Q M Q , r in /zT^Mpc and T in KeV. 




Figure 18. 2-D and 1-D marginalized posterior distributions for Mt(>2oo) 
(in h~^M Q x 10 14 ) and source flux densities (in Jys) within 5' from the 
cluster X-ray centroid for Abell 2111. 




Figure 20. 2-D and 1-D marginalized posterior distributions for Mx(r2oo) 
(in h~l Q x 10 14 M o ) and source flux densities (in Jys; see Tab. [5) within 5' 
from the cluster X-ray centroid for Abell 2146. 



expect this good agreement between these sets of results. We have 
further investigated the effect of sources in this cluster and have 
found a slight degeneracy between the cluster mass and the flux 
density of the source lying closest to the cluster centre - see Fig. 
[20lF - which had not been seen for the brighter, » 6mJybeam _1 
source lying a few arcseconds a way from the cluster centroid. 

Chandra data analysed bv lRussell et alJ (1201 Oh have revealed 
that Abell 2146 is undergoing a rare mer ger event similar to that 
of "Bullet-cluster" (IMarkevitch et al.ll2002r) . with two shock fronts 
with Mach numbers M « 2, and strong non-uniformities in the 



temperature profile. Note the different - essentially 90° - orien- 
taions between the X-ray and the SZ extensions. To understand this 
we have to consider collision geometry, mass ratio and, especially, 
time of snapshot since the merger start - see Sec. |8] 

6.15 Abell 2218 

Results for Abell 2218 are shown in Fig. [21] There is substantial 
radio emission towards Abell 2218, most of which is subtracted 
from our maps to leave a 470 //Jy beam -1 positive feature to the 
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Figure 16. Results for Abell 2009. Panels A and B show the SA map before and after source-subtraction, respectively; a 0.6 kA taper has been applied to B. 
The box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] The smoothed Chandra X-ray map overlaid with contours from 
B is presented in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F shows 
the 1 and 2-D marginalized posterior distributions for source flux densities (in Jys) given in Tab.|5]and M T (r 20 o) (in h'^ x 10 14 M o ). In panel D M T is given 
in units of h^ Q x 10 14 M o and f g in h^ Q ; both parameters are estimated within r2oo- In E M g is in units of h^ Q M Q , r in /z^Mpc and T in KeV 
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Figure 17. Results for Abell 2111. Panels A and B show the SA map before and after source-subtraction, respectively; a 0.6 taper has been applied to B. 
The box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] The smoothed Chandra X-ray map overlaid with contours from 
B is presented in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F shows 
the higher-resolution source-subtracted map (no taper). In panel D Mj is given in units of h^ Q x 10 14 and f g in h~^ Q ; both parameters are estimated within 
r2oo- In E M g is in units of h~Q Q M Q , r in /z^Mpc and T in KeV 
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Figure 19. Results for Abell 2146. Panels A and B show the SA map 
before and after source subtraction, respectively; a 0.6 k/l taper has been 
applied to B. The box in panels A and B indicates the cluster SZ cen- 
troid, for the other symbols see Tab. [4] The smoothed Chandra X-ray 
map overlaid with contours from B is presented in image C. Panels D 
and E show the marginalized posterior distributions for the cluster sam- 
pling and derived parameters, respectively. In panel D Mj is given in 
units of h^Q Q x 10 14 M© and f g in h^ Q ; both parameters are estimated 
within r2oo- in E M g is in units of h^Q Q M Q , r in /zjJ Mpc and T in KeV. 
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West of the decrement, which could be extended emission. Rudnick 
et al. detect diffuse emission from a radio halo with a flux of 0.05 Jy 
within a 500 kpc radius at 327 MHz, from which one might expect 

ign_ 
e.g jHanischUl98Gh . 

Several observations in the X-ray (e.g ., IMarkevitchl 



a < 200 ijl Jy s i gnal a t 16 GHz (for a typical halo spectral index, see 



1997: iGovoni et all |2004 iMachacek etall 12002), optical (e.g., 



Girad 



& Mez zetti 2001), SZ (e.g.. I Jones et al and lensing 

n 



Squires etal.lll996l and ISmith et all 120041) have suggested 



(e.g., 

that Abel l 2218 is a comp lex, disturbed syst em. High-resolution 
ROSAT (IMarkevitchl Il997h and Chandra jGovoni et all |2004 
IMachacek et alj|2002h data show signs of sub structure, particularly 
on small scales. Moreo ver, lensing studies by Is quires et all (119961) 
and lSmith et al.l j2004l) have revealed a bi-modal mass distribution 
and associated elongated structures in the mass distribution. 
Abell 2218 also sh ow s signs of strong temperature variations 
dGovoni et al.ll2004l and lPratt et ai1l2004l) . All of these results are 
indicative that the cluster is not relaxed. 

SZ observation s towards Abell 2218 have been made with 
the Ryle Telescope ( Jones et al . 2005) at 15 GHz, at 36 GHz using 
the Nobey ama Telescope ([Tsuboi etal.l ll998) and with OCRA-p 
at 30 GHz (lLancastej|2007h. Earlier SZ observations towards this 



cluster include iBirkinshaw et all Jl9 81). Birkinsh aw et all jl984l) . 



19911) . Uones et all (119931) and 



[Partridge et all dl987l) . iKlein et al 
Birkinshaw! (Il994l) . 

Pratt et al. find from XMM-Newton data that T(r) fal l s from 
8keV near the centre to 6.6 keV at 700 kpc. IZhang et al.l ([2008) 
calculate a cluster mass estimate from the XMM-Newton data; us- 
ing h 10 = 1.0, they obtain M T (r 500 ) = 4.2 + 1.3 x 10 14 M o and 
/ g (r 500 ) = 0.15 ± 0.09. We find M T (r 50 o) = 2.7 ± 0.6 X l0 u h^ Q M Q . 

The Chandra X-ray image shown in Fig. [21] C appears to be 
extended along the N-S direction on arcminute scales and along 
the ^SE-NW direction on scales « 2'. On the other hand, the dis- 
tribution of the X-ray signal on larger scales, « 3', tends to be more 
circular. On the untapered, source- subtracted SA map, Fig. |2T1F. 
the SZ signal towards Abell 2218 is clearly extended. There is no 
significant degeneracy between the cluster mass and the source flux 
densities, as seen from Fig. [22] 



6.16 Abell 2409 

We detect a 12cr 5A SZ effect towards Abell 2409 in the tapered, 
source- subtracted SA maps, Fig. [23] B. Despite the high SNR we 
are not able to obtain sensible parameter estimates for this cluster. 
As shown in Fig. [23] A, the effect of some emission close to the 
pointing centre is to give the decrement a shape that cannot be well 
approximated by a spherical /^-profile with free shape parameters. 
The parameter estimates from McAdam are thus not reliable and 
we present only the AMI SA map. Fixing the shape of the profile 
can improve the fit to this cluster. Cluster parameters for Abell 2409 
from AMI data have been obtained using a gNFW parameterization 
- see lPlanck and AMI Collaborations et all ( 1201 j) . 

The nature of the residual emission around the cluster is uncer- 
tain. Pointed LA observations towards the location of these sources 
of positive flux were made in an attempt to detect possible sources 
lying just below our detection threshold. Despite the noise at these 
locations on the LA map reaching « 50// Jy beam -1 , no additional 
sources were detected; it seems likely that this is (at this resolution) 
extended emission with relatively low surface brightness. However, 
no evidence for extended emission was found in either the NVSS 
1.4 GHz or in the VLSS 74 MHz maps. 




M T ( r 200) 



Figure 22. 2-D and 1-D marginalized posterior distributions for Mx(r2oo) 
(in h~l Q x 10 14 M o ) and source flux densities (in Jys) within 5' from the 
cluster X-ray centroid for Abell 2218 (Tab.[5}. 



6.17 RXJ0142+2131 

The maps and parameters for RXJ0142+2131 are presented in 
Fig-Ell The source environment is not expected to contaminate 
our SZ detection, with the brightest source having a flux density 
of « 2 mJy and lying several arcminutes away from the cluster cen- 
troid; residual emission after source subtraction is seen on the SA 
maps at the lcr level. The composite image of the SZ and X-ray 
data reveals good agreement between the emission peaks of these 
two datasets. 

A photometric a nd spectroscopic study of RXJ0142+2131 
bv lBaars et al ] d2005h finds the velocity dispersion of this cluster 
(<r x = 1278 + 134kms _1 ) to be surprisingly large, given its X-ray 
luminosity (Tab. [I}. This study indicates that galaxies in this cluster 
have older luminosity-weighted mean ages than expected, which 
could be explained by a short increase in the star formation rate, 
possibly due to a cluster-cluster merger. Moreover, RXJ0142+2131 
shows signs of not being fully virialized since the brightest cluster 
galaxy was foun d to be displaced by 1000 km s _1 from the systemic 
cluster redshift. lOkabe et all d2010h fitted an NFW profile for the 
mass density to Subaru/Suprime-Cam data and assumed a spherical 
geometry for the cluster to derive M T (r 500 ) = 2.85+fj g x l0 M h^M Q 
and M T (r 20 o) = 3.86+g;|| x lO^fc^M© (using h 72 = 1.0). From 
our analysis, we find M T (r 50 o) = 1.7 + 0.6 x 10 14 



/i loo M and 



M T (r 2 oo) : 



3.7^ x lO^ooM©. 



6.18 RXJ1720.1+2638 

Results for RXJ1720. 1+2638 are given in Fig. [25] At 16 GHz the 
source environment around the cluster is challenging: in our LA 
data we detect a 3.9 mJy source at the same position as the cluster, 
and several other sources with comparable flux densities within 4' 
from the cluster centre. The difficulty of modelling this system is 
clear from the degeneracies between some of the source flux densi- 
ties and the cluster mass (Fig.[25]F). However, we always recover 
a similarly asymmetric SZ decrement. 
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Figure 21. Results for Abell 2218. Panels A and B show the SA map before and after source subtraction, respectively; a 0.6 kA taper has been applied to B. 
The box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] The smoothed Chandra X-ray map overlaid with contours from 
B is presented in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. In panel 
D Mt is given in units of h~Q Q x 10 14 M© and f g in h~Q Q ; both parameters are estimated within r2oo- In E M g is in units of h~Q Q M Q , r in /z^Mpc and T in KeV. 
F shows the higher-resolution source-subtracted map (no taper). Contours of the map in panel C are not the same as in panel B.; they range from -1.388 to 
-0.188 mJybeam" 1 in steps of +0.15 mJybeam" 1 . © 0000 RAS, MNRAS 000, 000-000 
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Abell 2409 



A 




Figure 23. Results for Abell 2409. A: Source- subtracted SA map produced using a 0.6-k/l taper. B: SA map from A. overlaid onto Chandra X-ray image. 
Contours increase linearly in units of osa • 



RXJ1720. 1+2638 has been stud ied by iMazzotta et aD (1200 ll) 
and [Mazzotta & Giac intuccl (2008) through Chandra observa- 
tions. This cool-core cluster has two cold fronts within 100" 
of the X-ray centroid and a regular morphology away from 
the core region; the authors attribute the dynamics of this 
clust er to the sloshing scen a rio, i n agreement with later work 
by lowers. Nulsen. & Couchl (1201 lh using optical spectroscopy. 
Merger activity has also been suggested by Okabe et al. whose 
weak lensing data reveal a second mass concentration to the North 
of the main cluster, while the analysis of SDSS data bv lMiller et al.1 
J2005I) finds no evidence of substructure. Our data reveal a strong 
abundance of radio emission towards this cluster , including some 
extended emission, which might support the suggestion in Maz- 
zotta et al. (2001) that this cluster contains a low-frequency radio 
halo that did not disappear after the merger event. 

Mazzotta et al. 2001 determined the mass profile for the clus- 
ter assuming hydrostatic equilibrium to be M T (< r - lOOOkpc) = 
5^ x lO^h^M©. We find M T (r 50 o) = 1.2 + 0.2 x l0 u h^ M o . 



6.19 Zw0857.9+2107 

We report a null detection of an SZ signal towards this clus- 
ter, despite the low noise levels on our SA maps and a seem- 
ingly benign source environment. We reached a noise level (lcr) of 
97/Jy beam -1 on the LA map (Fig. [26] middle panel) and found no 
evidence for sources below our 4<r LA detection threshold. We detect 
a 1.4 mJy radio source at the location of the peak X-ray signal (see 
the electronic version in the ACCEPT Chandra data archive for a 
higher resolution X-ray image) but we seem to be able to subtract 
it well from the SA maps. 

Zw0857. 9+2107 is not a well-studied cluster. There are two 
temperature measurements for the cluster ga s in Zw0857.9+210 7 
from the ACCEPT Chandra data archive dCavagnolo et al.ll2000h : 
T * 3 + 4keV between * 10 < r < 100 kpc and T * 4.2 + 2.2 keV 
between « 100 < r < 600 kpc. One might expect the average tem- 



perature for the cluster to be even lower at larger radii, such that 
T(r2oo) < 3 keV. The absence of an SZ signal could be explained 
by a sharp radial drop in T or, perhaps, this cluster is particularly 
dense and compact such that it is X-ray bright but does not produce 
a strong SZ signal on the scales AMI is sensitive to (Alastair Edge, 
private communication). Fig. [27] illustrates what the marginalized 
parameter distributions look like for non-detections such as this. 



6.20 Zwl454.8+2233 

We detect no SZ effect in the AMI data towards Zwl454.8+2233, 
despite the low noise levels of our SA maps. We detect several 
sources close to the cluster centre, including ones with a flux den- 
sity of 1.64 mJy, 1.55 mJy and 8.4 mJy (at 13", 1.8' and 4.3' away 
from the pointing centre, respectively). The SA maps and derived 
parameters are shown in Fig. [27] The derived parameters for this 
non-detection are as expected: we find that M T (r 2 oo) approaches 
our lower prior limit and that M g shows similiar behaviour (see 
e.g., Fig. [27]). 

Zhang et al. found M T (r 500 ) = 2.4 + 0.7 x 1O 14 M using XMM- 
Newton, assuming isothermality, spherical symmetr y and /z 70 = 1. 
Chandra X-ray observations by Bauer et al. (2005) suggest the 
cluster has a cooling flow and lVenturi et all d2008) find from 610- 
MHz GMRT observations that the cluster has a core-halo radio 
source. 



7 SOURCE-SUBTRACTION SIMULATION 

Extracting robust cluster parameters for a system like Abell 2146 
with bright sources lying at or very close to the cluster is extremely 
challenging. Many factors can affect the reliability of the detec- 
tion and of the recovered parameters. Aside from model assump- 
tions, other important factors are: the SNR of the decrement in 
our maps, the wv-coverage, the size of the cluster and the distance 
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Figure 24. Results for RXJ0142+2131. Panels A and B show the SA map before and after source subtraction, respectively; a 0.6 kA taper has been applied to 
B. The box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab.|4] The smoothed Chandra X-ray map overlaid with contours 
from B is shown in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F 
shows the 1 and 2-D marginalized posterior distributions for source flux densities (in Jys) within 5' of the cluster SZ centroid (see Tab. [5} and Mjfooo) (in 
h^ x 10 14 M o ). In panel D Mj is given in units of h~Q Q x 10 14 and f g in h^ Q ; both parameters are estimated within r2oo- In E M g is in units of h^ Q M , r in 
/z7i n Mpc and T in KeV. 
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Figure 25. Results for RXJ1720. 1+2638. Panels A and B show the SA map before and after source subtraction, respectively; a 0.6 kA taper has been applied 
to B. The box in panels A and B indicates the cluster SZ centroid, for the other symbols see Tab. [4] The smoothed Chandra X-ray map overlaid with contours 
from B is shown in image C. Panels D and E show the marginalized posterior distributions for the cluster sampling and derived parameters, respectively. F 
shows the 1 and 2-D marginalized posterior distributions for source flux densities (in Jys) given in Tab.|5]and Afrfooo) (in h^ x 10 14 M o ). In panel D is 
given in units of h^ Q x 10 14 M o and f g in h^ Q ; both parameters are estimated within r2oo- in E M g is in units of h^ Q M Q , r in /z^Mpc and T in KeV. 
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Zwl454.8+2233 
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Figure 27. The null detection of Zwl454. 8+2233 in SZ. Panel A shows the SA map before subtraction, which reveals the challenging source environment 
towards this cluster. The SA map after source subtraction is shown in panel B; no convincing SZ decrement is visible. Image C shows the Chandra X-ray map 
overlaid with SA contours from panel B. Panels D and E show the distributions for the sampling and derived parameters respectively; such distributions are 
consistent with a null detection. In panel D Mj is given in units of h^ x 10 14 M o and f g in h^ ; both parameters are estimated within r2oo- In E M g is in units 
of h~Q M Q , r in /z^Mpc and T in KeV. 
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Table 9. Source parameters for the three simulated sources. 



9 h 01 m 10 s 01 m 00 s 00 m 50 s 00 m 40 s 00 m 30 s 

Figure 26. Top panel: Source- subtracted SA map for Zw0857. 9+2107 pro- 
duced using a 0.6-kA taper. The contours increase linearly in units of ctsa- 
Middle panel: LA signal-to-noise map. Contours start at 3cr and increase 
linearly to 10cr, where <x = 97rawJy Beam -1 . Bottom panel: SA contours 
overlaid onto the Chandra X-ray image. The SA contours are the same as in 
the upper panel. 
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Figure 28. SA contour map of simulated data containing thermal noise + 
confusion noise + CMB + cluster + resolved point sources. Contours in- 
crease linearly in units of <tsa- 



of the sources from the cluster, their flux-densities and their mor- 
phologies. From Sec. [6] one can appreciate that at 16 GHz the SZ 
signal is potentially strongly contaminated by radio sources. We 
have examined some of the effects of these sources in a controlled 
environment through simulations. For this purpose, we generated 
mock visibilities between hour angles -4.0 to 4.0, with an RMS 
noise per channel per baseline per second of 0.54 Jy. Noise con- 
tributions from a CMB realisation and from confusion from faint 
sources lying below our subtraction limit were included; for the 
former we used a ACDM model and for the latter we integrated the 
IOC LA source counts from 10// Jy to 300yt/Jy. A cluster at z = 0.23 
was simulated using an isothermal ^-profile to model the gas dis- 
tribution, with a central electron density of 9 x 10 3 m -3 , p = 1.85, 
r c = 440/z~ 1 kpc and T = 4.8 keV. Integrating the density profile 
out to r 2 oo (Eq. HJ assuming a spherical cluster geometry yields 
M g (r 2 oo) = 6.25 x 10 13 /z;^M o . From the virial M-T relation given 
in Eq.[4]M T (r2oo) = 5.70x lO^^M© and using these two estimates 
and Eq.|5] we find / g (r 2 oo) = 0.11/z^J. 

Three point sources were included into the simulation. Their 
positions, flux-densities and spectral indices are given in Tab. [9] 
The map of the data is shown in Fig. [28] 

The data for the simulation were run through the same analysis 
as described in Sec. [5] In this case the source priors were centred on 
the simulated values (Tab. [9} and the cluster priors were the same 
of those in Tab. [3] with the delta-prior on z set to 0.23. The 1-D and 
2-D marginalized posterior distributions for the sampling parame- 
ters are presented in Fig. [29] It can be seen that the cluster position 
and gas fraction are recovered well by the sampler; the core radius 
and p cannot be constrained by AMI data alone, thus, as expected, 
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Table 10. Median, mean and standard deviation for Msz,sph/Msz,eiiip 
within r2oo and r^oo and TAMi,sph/TAMi,eMp- Data for all clusters in Tab. 
[15] were included, except for Abell 1758a and b. Ratios for each cluster at 
these two overdensities are given in Tab. [15] 

Median Mean Standard deviation 

Msz,sph/M S z,eiiip within r 500 0.96 0.96 0.16 

Msz, S ph/M S z,eiiip within r 200 0.97 0.99 0.16 

TAMI,sph/TAMI,ellip 0.98 0.98 0-10 



the agreement between the input and output mean values for these 
paramaters is poor; the total cluster mass, on the other hand, is very 
well-constrained and the recovered value is consistent with the in- 
put value. Hence, despite the challenging source environment, and 
the degeneracies between the cluster mass and the source flux den- 
sities, our analysis is able to provide robust cluster mass estimates. 



8 DISCUSSION 

Of the 20 target clusters, we have detected SZ towards 17, all 
of which are resolved, and with "peak" detections between 5cr SA 
and 23ct S a- The analysis has produced robust parameter extrac- 
tion for 16 of the 17 - this was not possible for Abell 2409 be- 
cause of nearby extended radio emission that distorts the SZ signal 
and gives an unacceptable fit for a spherical /?-model. The three 
null detections are of Abell 1704 (difficult source environment), 
Zw0857. 9+2107 (it is unclear to us why we have not detected this), 
and Zw 145 8. 8 +223 3 (difficult source environment). 

8.1 Cluster morphology and dynamics 

The images frequently show significant differences in position of 
the SZ peak (and of the SZ centroid) and the X-ray peak, indicating 
that the densest part of a cluster is not at the centre of the large-scale 
gas distribution. In Abell 773 and Abell 2146, both mergers, there 
is evidence of SZ extension perpendicular to the X-ray emission. 
Abell 1758a and Abell A1758b are both major mergers and there is 
a hint of an SZ signal between a & b. Unlike what one might naively 
expect, there are cases of SZ extensions in non-mergers and cases 
of near-circular SZ map structures in mergers. 

To attempt to quantify the cluster morphology from the AMI 
data, we ran our analysis with an ellipsoidal model for the cluster 
geometry. This model simply fits for two additional parameters: an 
ellipticity parameter, r/, which is the ratio between the semi major 
and semi minor axes and an angle 6 measured anticlockwise from 
the West; these val ues are given in Tab. ITT1 For further de t ails o n 
this model see e.g. J AMI Consortium: Hurley-Walker et"aD d201 lh . 

As a check that switching from spherical to ellipsoidal SZ 
analysis does not itself introduce significant bias in mass, we give 
in Tab. [TOl the ratios M S z,sphlM S z,eiiip within r 20 o and r 500 and 
TAMi,sph/TAMi,eiii P - no significant bias is evident. Of course, else- 
where in this paper we use spherical SZ estimates because the X- 
ray and almost all the optical total cluster mass estimates also as- 
sume spherical symmetry. 

Tab. [TT] also includes other possible indicators of dynamical 
state. The presence of cooling cores (CC) is associated with re- 
laxed clusters since it is widel y accepted that merger e vents tend 
to disrupt cooling flows, e.g. jMcGlvnn & Fabianl (11984). We have 
used Chandra data from the ACCEPT dat abase, where availab le, to 
compute three CC indicators described in lHudson et al.ld2010h : the 
central entropy, the central cooling time and the ratio of (approx- 
imately) the central cluster temperature to the virial temperature; 



Tab. [TT] also includes other assessments of dynamical state that we 
have found in the literature. 

The projected separation of the brightest cluster galaxy (BCG) 
and the peak of the X-ray emission has been shown to corre- 
late with the dynamic al eq uilibrium state of the the host cluster 
(iKatavama etaJ 120031 and ISanderson et all [2009h . Similarly, the 
offset between the SZ centroid and the X-ray peak can also be 
a diagnostic for cluster disturbance. For this purpose, the separa- 
tion between the AMI SZ centroid, X-ray peak cluster position and 
the position of the BCG are given in Tab. Q21 in Tab. [13] some 
sample statistics are provided. Large offsets betw een these mea- 
surements have been reported in observations (e.g.. | Massardi et al] 
l201(j|Menanteau et alJ[201lL andlKorngut et alfcoilh and in sim- 
ulations (e.g., dMolnar. Hearn. & Stadelll2012h ). 

Examination of Tabs. ITTIIT31 indicates that even for well- 
studied clusters there are conflicting indications as to whether the 
cluster is a merger or not, e.g., Abell 773 does not appear to have a 
CC, has high degree of ellipticity, the X-ray and SZ signals appear 
to be oriented quasi-perpendicularly to each other and yet the rela- 
tively small position offsets in Tab. [12] might suggest the cluster is 
relaxed. 

8.2 SZ temperature, large-radius X-ray temperature, and 
dynamics 

In Fig.[30lwe compare the AMI SA observed cluster temperatures 
within r 2 oo (?ami) with large-radius X-ray values (T x ) from Chan- 
dra or Suzaku that we have been able to find in the literature. We 
use large-radius (« 500 kpc) X-ray temperature values to be con- 
sistent with the angular scales measured by AMI. (For Abell 611 
we have plotted two X-ray values from Cha ndra data - one from 
the ACCEPT archive dCavagnolo et alj200QL which is higher than 
our AMI SA measurement, and a second X-ray measurement from 
Chandra (Donnarumma et al.), which is consistent with our mea- 
surement). 

There is reasonable correspondence between SZ and X-ray 
temperatures at lower X-ray luminosity, with excess (over SZ) X- 
ray temperatures at higher X-ray luminosity. The mean, median 
and standard deviation for the ratio of T AM i/T x were found to 
be 0.7, 0.8 and 0.2, respectively, when considering all the clus- 
ter in Fig. [30] )except for Abell 1758a, due to it being a com- 
plex double-merger). The numbers are obviously small, but the 
two systems that are strong mergers by clear historical consensus - 
Abell 773 and Abell 1758a - are unambiguously clear outliers with 
much high er large-radius X-ray temperatures than SZ temperatures. 
ISmith et al] d2004h investigate the scatter between lensing masses 
within < 500 kpc with Chandra X-ray temperatures averaged over 
0. 1 - 2 Mpc for ten clusters and also find that disturbe d systems have 
higher temperatures. However, Marr one et al measure the 

relationship between SZ-7 sph and lensing masses within 350 kpc for 
14 cluster s and find no segregation between d isturbed and relaxed 
systems. iKravtsov. Vikhlinin. & Nagail |2006) analysed a cluster 
sample extracted from cosmological simulations and noticed that 
X-ray temperatures of disturbed clusters were biased high, while 
the X-ray analogue of SZ-7 sph , did not depend strongly on cluster 
structure. Taken together, these results suggest that, even at small 
distances from the core, SZ-based mass (or temperature) is a less 
sensitive indicator of disturbance disturbance than is X-ray-based 
mass. 

Major mergers in our sample have large-radius X-ray temper- 
atures (at « 500 kpc) higher than the SZ temperatures (averaged 
over the whole cluster). This suggests that the mergers affect the 
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Table 11. Dynamical indicators: 6, the angle measured anticlockwise from the West, 77 the ratio between the semi major and semi minor axes (these values 
arise from fitting the SZ data with an elliptical geometry (see text)). Cooling core information: CC denotes the presence of a cooling core and NCC the lack of 
('-' means this information is not clear or not known); Corei is a result from this study obtained by using three CC indicators described in Hudson et al. 2010 
- the central entropy, the central cooling time and the ratio of approximately the central cluster temperature to the virial temperature, 7o/r vir , where all the 
data have been taken from the Chandra ACCEPT database; Core2 is cooling core information on the cluster available from other studies, f : the core type of 
Abell 1423 is unclear; the ratio of Tq/T^ taken from ACCEPT suggests it is not a cool-core cluster but, a CC cannot be ruled out due to the large uncertainty 
in the X-ray temperature measurements; the central entropy and cooling time are unclear. 



Cluster Name 


e 


1 


Corei 


Core2 


Abell 586 


136 ± 33 


0.73 ± 0.13 


NCC 


CC Allen (2000), NCC Marrone et al. (201 1) 


Abell 611 


79 ± 39 


0.80 ±0.12 


NCC 


NCC Marrone et al. (2011) 


Abell 621 


64 ± 61 


0.73 ± 0.13 






Aoeii / ID 


41 ± 1Z 


U.jV ± U.1U 




Allen (ZUUUj 


Abell 781 


132 ± 32 


0.70 ±0.13 






Abell 990 


109 ± 40 


0.78 ±0.13 






Abell 1413 


101 ±21 


0.75 ±0.12 


CC 


CC Allen (2000). Richard et al. (2010) 


Abell 1423(f) 


66 ±28 


0.70 ±0.14 


NCC? 


& CC? Sanderson et al. (2009) 


Abell 1704 








CC Allen (2000) 


Abell 1758a 


72 ±31 


0.73 ±0.14 


NCC 




Abell 1758b 


85 ±49 


0.77 ±0.13 






Abell 2009 


88 ±49 


0.78 ±0.12 






Abell 2111 


90 ±29 


0.77 ±0.12 


NCC 




Abell 2146 


126 ±4 


0.56 ± 0.05 




'Bullet-like merger' Russell et al. (2010) 


Abell 2218 


107 ± 80 


0.87 ± 0.07 


NCC 


NCC Richard et al. (2010) 


Abell 2409 








RXJ0142+2131 


87 ±41 


0.77 ±0.13 






RXJ1720. 1+263 


26 ± 12 


0.58 ± 0.07 


CC 


CC Richard et al. (2010) 


Zw0857.9+2107 










Zwl454.8+2233 








CC Bauer et al. (2005) 



Table 12. X-ray cluster position; SZ centroids from our analysis; position of the BCG from SDSS (the BCG was identified as the brightest galaxy in the central 
few hundred kpc from the cluster X-ray position. For clusters labelled with (*) the BCG could not be identified unambiguously. Entries filled with a '-' indicate 
there is no available information. 



Cluster Name 


BCG 




X-ray 




SZ 






Position offsets (") 






RA (Deg) 


Dec (Deg) 


RA (Deg) 


Dec (Deg) 


RA (Deg) 


Dec (Deg) 


SZ-X-ray 


SZ-SDSS 


X-ray-SDSS 


Abell 586 


113.0844 


31.6334 


113.0833 


31.6328 


113.0833 


31.6264 


23.0 


25.7 


4.5 


Abell 611 


120.2367 


36.0563 


120.2458 


36.0503 


120.7958 


36.0531 


10.1 


34.9 


39.4 


Abell 621 






122.8000 


70.0408 


122.7875 


70.0458 


48.5 






Abell 773 


139.4724 


51.7270 


139.4666 


51.7319 


139.4667 


51.7331 


4.0 


29.9 


27.3 


Abell 78 1(*) 


140.1073 


30.4941 


140.1083 


30.5147 


140.1000 


30.5314 


67.0 


136.7 


74.2 


Abell 990 


155.9161 


49.1438 


155.9208 


49.1439 


155.9125 


49.1369 


39.0 


27.9 


16.9 


Abell 1413 


178.8250 


23.4050 


178.8250 


23.4078 


178.8250 


23.3894 


66.0 


55.8 


10.2 


Abell 1423 


179.3222 


33.6110 


179.3416 


33.6319 


179.3375 


33.6189 


49.2 


62.2 


103.0 


Abell 1704 


198.6025 


64.5753 


198.5917 


64.5750 










39.0 


Abell 1758a(*) 


203.1189 


50.4697 


203.1500 


50.4806 


179.3375 


50.5264 


209.3 


209.3 


118.6 


Abell 1758b 










203.1250 


50.4003 


209.3 


209.3 




Abell 2009 


225.0833 


21.3678 


225.0811 


21.3692 


225.0875 


21.3553 


55.2 


47.5 


9.5 


Abell 2111 


234.9333 


34.4156 


234.9187 


34.4240 


234.9125 


34.4331 


39.4 


97.9 


60.8 


Abell 2146 






239.0291 


66.3597 


239.0250 


66.3589 


15.1 






Abell 2218 






248.9666 


66.2139 


248.9375 


66.2186 


106.1 






Abell 2409 


330.2189 


20.9683 


330.2208 


20.9606 










28.5 


RXJ0142+2131 






25.51250 


21.5219 


25.51667 


21.5303 


33.6 






RXJ1720. 1+2638 


260.0418 


26.6256 


260.0416 


26.6250 


260.0333 


26.6125 


54.0 


56.0 


2.1 


Zw0857.9+2107 


135.1536 


20.8943 


135.1583 


20.9158 










79.5 


Zwl454.8+2233 


224.3130 


22.3428 


224.3125 


22.3417 










4.3 



Table 13. Mean, standard deviation and median for the differences in X-ray and SZ cluster centroids and the position of the BCG from SDSS maps. Abell 
1758 (a and b) has been excluded from this analysis due to its exceptionally disturbed state. 





Mean (") 


Standard Deviation (") 


Median (") 


SZ - Xray 


43.9 


26.8 


43.6 


SZ - SDSS 


51.6 


35.3 


43.6 


Xray - SDSS 


27.9 


32.2 


35.7 
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Figure 29. Left: One and two-dimensional marginalized posterior distributions for the cluster sampling parameters from our simulation. Mj and f g are 
estimated within r 2 oo and Mj is given in units of xlO 14 . The green crosses in the 2-D marginals denote the mean of the distribution. Right: One and two- 
dimensional marginalized posterior distributions for the source flux densities and Mx(r 2 oo) for our simulation. Red lines indicate the mean of the marginalized 
distribution and the blue lines represent the input value. 



n 2 -weighted X-ray temperatures more than the n- weighted SZ tem- 
peratures and do so out to large radius. This is evidence for shock- 
ing or clumping or both at large radius in mergers. Indications that 
clumping at large r might h ave a significant impact o n X-ray re- 
sults have been found by e.g jKawaharada etall d2010h . who find a 
flattening of the entropy profile around the virial radius, c ontrary to 
the theoretical predictions (e.g.. | Voit Kav. & Brvanll2005l) . Hydro- 
dynamical simulations bv lNagai & Laul fcOllh have shown that gas 
clumping can indeed introduce a large bias in large- r X-ray mea- 
surements and could help explain the results by e.g., Kawah arada et 
al. It should be noted, however, that lMazzotta et alj (120041) expect 
X-ray temperatures to be lower than mass-weighted temperatures 
for clusters with temperature structure since the detectors of Chan- 
dra and XMM-Newton are more efficient on the soft bands, which 
leads to an up weighting of the cold gas. However, in simulations by 
iRasia et all 020121) mass-weighted temperatures were shown to be 
larger than X-ray temperatures for the vast majority of their clus- 
ters, particularly for very the most disturbed clusters in their sam- 
ple. 

Examination of Fig. [30l given Tab.[TT]is suggestive of another 
relation, again with obviously small numbers. Fig. [31] shows AMI 
cluster temperature versus large-scale X-ray temperature but with 
each cluster X-ray luminosity replaced by AMI ellipicity, 77, and 
its error (note that we have removed the two Abell 611 points be- 
cause of their apparently discrepant X-ray values). With one excep- 
tion (Abell 2146), the clusters with large-radius X-ray temperature 
^ 6keV have 77 values < 0.70, whereas the first two outliers to the 
right (RXJ1720. 1+2638 and Abell 773) have significantly smaller 
values of AMI ellipticity. The rightmost outlier (Abell 1758a) it- 
self has the ellipticity value 0.73 + 0.14 but this will be mislead- 
ingly high if we should instead be considering the ellipticity of the 
Abell 1758a+b taken as a merging pair. The true relationship be- 
tween SZ ellipticity and merger state is bound to be influenced by 
the c ollision geometry, the time since the start of the merger (Fig. 
1 in iNelson et alJuOllI illustrates how SZ 77 and 6 can vary with 



merger evolution), the mass ratio, and so on. Far more data, includ- 
ing data on clusters not selected by X-ray luminosity, are essential. 

8.3 Comparison of masses within r 50 o and within the virial 
radius (« r 20 o) 

The classical virial radius, « r 20 o, found is typically 1.2+0.1 Mpc. 
Values for M T (r 20 o) range from 2.0^ x 10 14 ^ 1 M o to 6.1 + 0.9 x 
l0 u h~^M Q and are typically 2.0-2.5X larger than M T (r 500 ). In Fig. 
[32l and [33l AMI mass estimates at two overdensity radii are com- 
pared with other published mass estimates. The scarcity of mass 
measurements at large r is apparent from these figures. 

• For M T (r 5 oo), there is good agreement between optical and 
AMI (HSE) mass estimates. In contrast, the X-ray (HSE) estimates 
tend to be higher, sometimes substantially so. 

• For M T (r 20 o), there is very good agreement between optical es- 
timates, the Suzaku X-ray (HSE) estimate, and the AMI (M-T) es- 
timates. Good agreement b etween AMI and optical masses has pre- 
viousl y been reported by lAMI Consortium: Hurley-Walker et all 

boiih . 

• From our sample we cannot determine whether the disagree- 
ment of masses is a function of radius. 

• The discrepancy between the X-ray and AMI masses for 
Abell 1413 is reduced at r 20 o, with the X-ray mass being larger than 
the AMI mass by « 50% at r 50 o and smaller than the AMI mass by 
« 10% at r 200 . 

• The largest discrepancies between mass measurements in SZ, 
optical and X-ray correspond to the strongest mergers within our 
sample but the X-ray masses are always higher than our SZ masses, 
even for the few relaxed clusters in our sample. Given that the lens- 
ing masses agree well with our SZ estimates, this might be an indi- 
cation of a stronger bias in masses estimated from X-ray data than 
from SZ or lensing data, especially for disturbed systems. How- 
ever, most recent simulations and analyses indicate that X-ray HSE 
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Large-radius X-ray temperature (keV) 

Figure 30. The AMI mean temperature within r2oo versus the X-ray temperature. Eac h point is labelled with the cluster name and X-ray luminosity. Most 
of the X-ray measurements are large-radius temperatures from the ACCEPT archive (Cav agnolo et ail 2000) with 90% confidence bars. The radius of the 
measurements taken from the ACCEPT archive are 400-600 kpc for Abell 586, 300-700 kpc for Abell 611, 300-600 kpc for Abell 773, 450-700 kpc for 
Abell 1423, 500-1000 kpc for A bell 2111, 450-550 for Abe ll 2218 and for RXJ1720. 1+2638 r = 550-700 kpc. The Abell 611* temperature is the 450-750 kpc 
value with 68% confidence bars ([Donnarumma et al. 2010). The Abell 2146 temperature measurement is from Russell et a l. 2010 (with 68% con fidence bars). 
The Abell 141 3 X-ray temperature i s es timated from the 700 - 1200 kpc measurements made with the Suzaku satellite (Hoshino et al. 2010), this value is 
consistent with lVikhlinin et aljEooil and lSnowden et al.ll2008l The ACCEPT archive temperature for Abell 175 8 A is 16+7 keV at r= 475-550 kpc, and with 
SZ temperature 4.5±0.5, is off the right-hand edge of this plot. Abell 611 has been plotted using dashed blue lines to emphasize that this cluster has two 
X-ray-derived large- r temperatures. The black diagonal solid line is the 1:1 line. 



masses are underestimated with respect to lensing masse s (e.g., 
iNagai et al.ll2007llMeneghetti et all201cLlRasia et al.ll2012h . 

8.3.1 Related results from the literature 

To illustrate some of the issues in mass estimation, we bring to- 
gether some of the other results in the literature. 

• X-ray and weak-l ensing masses 

Observational studies by Mahdavi et al. ( 2008), Zh ang et all ( 2008) 
and lZhang et al] (120 101) find systematic differences between X-ray 
HSE-derived and weak-lensing masses, with the lensing masses 
typically exceeding the X-ray masses. Madhavi et al. report a strong 
radial dependence for this difference, with weak-lensing masses be- 
ing ^3% smaller within r 2 5oo but « 20% larger within r 50 o than the 
X-ray masses, yet find no correlation between the difference level 
and the presence of cool cores. Zhang et al. (2010) find that X-ray 
masses seem underestimated by « 10% for undisturbed systems 
and overestimated by « 6% for disturbed clusters within r 500 . For 
relaxed clusters, they find the discrepancy is reduced at larger over- 
densities. The underestimate of HSE X-ray masses with respect 



to lensing masses has been widely produ ced in simulations (e .g., 
INagai et al.ll200llMeneghetti et alj2O10h and lRasia et aDl2012l . 

In Tab. [14] we follow Mahdavi et al. to calculate a weighted best- 
fit ratio of two mass estimates at different overdensities for different 
data. The simulations by Rasia et al. and Meneghetti et al. yield 
signi ficantly lower M x /M WL at r 50 o than the observational data. 
Siiacki et all d2007h suggest that a higher incidence of temperature 
substructure in the simulations might be responsible for this effect. 
It is interesting to see how the mass agreement for the study by 
Zhang et al. seems to weaken when excluding disturbed systems. 
What is very different from the literature is that we find HSE X-ray 
masses to be consistently higher than our HSE SZ masses within 
r 50 o. Modelling our clusters with an elliptical model for the cluster 
geometry does not substantially improve the agreement. 

• SZ Y with X-ray and lensing masses 

iBonamente et al] {|2012h find good agreement between r sph (r 5 oo) es- 
timated from a joint S Z and X-ray analysis and from SZ d ata alone, 
in support of results by Planck Collaborati on et al .1 d201 lh . For their 
sample of massive, relaxed clusters there appears to be no signif- 
icant systematics affecting the ICM pressure measurements from 
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Figure 31. Plot analogous to Fig.[30lbut with the X-ray luminosity values replaced by SZ 77 (ellipticity). 



X-ray or SZ data. But, of course, this result might not be repro- 
d uced for a sample of d isturbed clusters. 

iMarrone et al 1 d2009h measure the scaling between 7 S z and weak 
lensing mass measurements within 350 kpc (« r 40 oo-800o) for 14 
LoCuSS clusters. They find it behaves consistently with the self- 
similar predictions, has considerably less scatter than the relation 
between lensing mass and T x and does not depend strongly on the 
dynamical state of the cluster. They suggest SZ parameters derived 
from observations near the cluster cores may be less sensitive to 
the complicate d physics of t hese r egions than those in X-ray. A 
later study by IMarrone et al.l d201lh comparing two F S z - M scal- 
ing relations using weak-lensing masses and X-ray (HSE) masses 
at r 2 5oo, nooo and r 50 o indicates the latter has more scatter and is 
more sensitive to cluster morphology, with the mass estimates of 
undisturbed clusters exceeding those of disturbed clusters at fixed 
^s P h by & 40% at large overdensities. However, this division is not 
predicted by comparing SZ and true masses from simulations and 
is could due to the use of a simple spherical lens model. Moreover, 
recently, iRasia et all d2012h have shown through simulations that 
selecting relaxed clusters for weak-lensing studies based on X-ray 
morphology is not optimal since there can be mass from, e.g., fila- 
ments not associated to X-ray counterparts biasing the lensing mass 
estimates even for systems which appear to be regular in X-rays. 

• Simulations 

Simulations of cluster mergers have shown these events generate 
turbulence, bulk flows and complex temperat ure structure, all o f 
which can result in cluster mass biases (e.g., IPoole etal.lE007h . 



Predominantly, simulation s indicate that X-ra y HSE masses tend 
to be underestimated (e.g. jKrause et al.|[2012h particularly in dis- 
turbed clusters, though the amount of the bias varies depending on 
the the simulation details, particularly on the physical processes 
taken into consideration. Projections effects, model assumptions 
and the dynamical state of the cluster are some of the factors affect- 
ing h ow well the true cluster mass can be m easured. As shown by 
e.g jTakizawa. Nagino. & Matsushital J2010h , even mass estimates 
for spherical X- ray systems are not always recovered well. Recent 
simulations by iNelson etaD fcOllh have investigated in detail the 
evolution of the non-thermal support bias as function of radius and 
of the merger stage. They reveal a very complex picture: the HSE 
bias appears to vary in amplitude and direction radially and as the 
merger evolves (and the shocks propagate through); for the most 
part, the HSE bias leads to an underestimate for the mass, there are 
times when it has the opposite effect. 

From simulations there appear to be two main, competing ef- 
fects that can lead to a mass bias from the effects of a merger. 
Firstly, the mer ger event can boost the X-ray luminosity and tem- 
perature (e.g., Rick er & SarazinlEoOll) such that if the cluster is 
observed during this period its X-ray mass will be overestimated. 
Secondly, the increase in non-thermal pressure support during the 
merger can lea d to X-ray (HSE) cluster m asses being underes- 
timated (e.g., lLau. Kravtsov. & Nag ai 2009). The cluster sample 
derive d from simulations studied bv lKravtsov. Vikhlinin. & Nagail 
(2006) showed that the X-ray temperatures were biased high for 
disturbed clusters, unlike Y x , the product of the gas mass and tem- 
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Table 14. Best-fit mass ratios calculated following Mahdavi et al. 2008. 
R12 are the results from simulations by Rasia et al. 2012, ME 10 are the 
simulations from Meneghetti et al. 2010, Z10 from Zhang et al. 2010 and 
MA10 from Mahdavi et al 2008. For our results we have used for simplicity 
sph to denote our SZ masses derived using a spherical geometry and ellip 
when assuming an elliptical model. We have excluded Abell 1758 (A and 
B) from the analysis, given its abnormally disturbed and complex nature. 

^500 ^200 



M X /M WL 



R12- full sample 


0.75 


± 


0.02 


R12- regular clusters 


0.75 


± 


0.04 


ME10- full sample 


0.88 


± 


0.02 


Z 10- full sample 


0.99 


± 


0.07 


Z10- relaxed 


0.91 


± 


0.06 


MA10- all 


0.78 


± 


0.09 


This work 









Mx/M S z,sph 1.7 ±0.2 

Mx/M S z,eiiip 1.6 ±0.3 

M SZ , S ph/M WL 1.2!jj| 1.0 ±0.1 

Msz,eiii P /M WL 1.2 + 0;jj 0.9 ±0.1 



perature as deduced from X-ray observations (the X-ray analogue 
of the SZ Y) which did not appear to depend strongly on cluster 
structure. 



9 CONCLUSIONS 

We observe 19 LoCuSS clusters with L x > 1 x 10 37 W (h 50 = 1.0) 
and present SZ images before and after source subtraction for 16 
of them (and for Abell 1758b, which was found within the FoV of 
Abell 1758a). We do not detect SZ effects towards Zwl458.8+2233 
and Abell 1704, due to difficult source environments, nor towards 
Zw0857.9+2107, for reasons unclear to us. We have produced 
marginalized posterior distributions at r 50 o and r 2 oo for 16 clusters 
(since Abell 2409 can not be fitted adequately by our model). 

• Measurements of M T (r 2 oo) are not common in the literature 
but are very important for testing large-radius scaling relations 
and understanding the physics in the outskirts of clusters. Conse- 
quently, the 16 measurements presented here, from a sample with 
narrow redshift-range, represent a significant increment to what al- 
ready exists. 

• For the clusters studied, we find values for M T (r 2 oo) span 2.0- 
6.1 + 0.9 x 10 14 h^M Q and are typically 2-2.5 times larger than 
M T (r 5 oo); we find r 2 oo is typically 1.1 + O.l/^ 1 Mpc. 

• AMI measurements of M T (r 50 o) are consistent with pub- 
lished optical results for 3 out of 4 clusters in our sample, with 
the weighted best-fit ratic|3 between AMI SZ masses and lensing 
masses being 1.2^ 3 within r 50 o and 1.0 + 0.1 within r 20 o- They are 
systematically lower than existing X-ray measurements of M T (r 50 o) 
for 6 clusters with available X-ray estimates and are only consis- 
tent with one of these measurements. The more discrepant masses 
correspond to the stronger mergers of the sample. The ratio of the 
X-ray masses to the AMI SZ masses is 1.7 + 0.2 for the sample. 
The agreement with optical measurements improves for M T (r 20 o), 
though there are few data. We have investigated the AMI vs X- 
ray discrepancy by comparing T AMI estimates with T x estimates, 
when available, at r « 500 kpc. There tends to be good agreement 

9 With the exception of Abell 1758a+b 



in less X-ray luminous clusters and in non-mergers but large -radius 
T x can be substantially larger than T AM i in mergers. This explains 
why some X-ray mass estimates are significantly higher than the 
AMI estimates: the use of a higher temperature will give a con- 
sequently higher mass in the hydrostatic equilibrium model used. 
Another implication of a higher large-radius T x than T A mi (given 
the respective n 2 and n emission weightings) is that, even at around 
r 50 o, the gas is clumped or shocked or both. There is a clear need 
for more large-scale measurements. 

• We have investigated the effects of our main contaminant, ra- 
dio sources, by searching for degeneracies in the posterior distri- 
butions of source flux densities for sources within 5' of the clus- 
ter SZ centroids. We find small or negligible degeneracies between 
source flux densities and cluster mass for all clusters, with the ex- 
ceptions of Abell 781, Abell 1758a and RXJ1720. 1+2638, which 
have sources with flux densities of 9, 7 and 7mJy at < 2' from 
the cluster SZ centroids. By simulating a cluster with a challeng- 
ing source environment, we have shown that our AMI analysis can 
approximately recover the true mass, even in a degenerate scenario. 

• We often find differences in the position of SZ and X-ray 
peaks, with an average offset of 35", a median of 34" and a sam- 
ple standard deviation of 24" for the entire sample (excluding 
Abell 1758a+b), confirming what has been seen in previous obser- 
vational studies and in simulations. We emphasize that our sample 
size is small, but we find no clear relation (except for Abell 1758) 
between position difference and merger activity. There is however 
an indication of a relation between merger activity and SZ elliptic- 
ity. 

• We have analysed the AMI data for two clusters: Abell 611 
and Abell 2111, with a J3 parameterization and with five gNFW 
parameterizations, in cluding the widely used lArnaud et al.l 02010) 
"universal" and the iNagai et all d2007h ones. This has revealed 
very different degeneracies in r sph (r 5 oo) - r 50 o for the two types of 
cluster parameterization. For both clusters, the ft parameterization, 
which allows the shape parameters to be fitted, yielded stronger 
constraints on r 50 o than any of the gNFW paramaterizations. The 
Nagai et al. and Arnaud et al. gNFW parameters produced consis- 
tent results, with the latter giving slightly better constraints. Setting 
the gNFW parameters to different, but reasonable, values altered 
the degeneracies significantly. This illustrates the risks of using a 
single set of fixed, averaged profile shape parameters to model all 
clusters. 
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Figure 32. Comparison of AMI Mj(r '500) measurements with others. Methods use d for estim ating Mj(rsoo) are giv en in the legend. The line of grad ient one 
has been included to aid the comparison. The references are as follows: Abell 586 (Okabe et al.ll2010h : Ab ell 611 {Okabe et al.ll2010h : Abe ll 773 (|Zhang et al.1 
|201Q|) : Abell 781 dSehgal et al.ll2008l andlZhang et all2010l) : Abell 1413 dZhang et all2010h . Abell 1758A dZhang et all2010h . Abell 2218 dZhang et al.l201Qh 
and RXJ0142+2131 dOkabe et al.ll201(J) . AMI values are given in Tab. [6] These were the M(>5oo) from X-ray and weak lensing data that we found in the 
literature. 
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Figure 33. Comparison of AMI Mt(>2oo) measurements with others. Methods used for estimating Mj(r 500) are given in the lege nd. Mass is give n 
in units of xlO 14 M . The lin e of gradient one has bee n included to aid comparison. The reference s are as follow s: Abell 586 dOkabe et al.l 2010); 
Abell 611 |bkabe et al.l I2010L iRomario et al.l l20ld and lAMI Consortium: Hurlev- Walker et all 1201 lh : Abell 1413 dHoshino et all bold) : Abell 2111 
dAMI Consortium: Hurley- Walker et al.ll201ll) and RXJ0142+2131 dOkabe et alj|2010h . AMI values are given in Tab.|g] 
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Table 15. Comparison of cluster masses at and r2oo for a spherical and an elliptical model for the cluster geometry. Ratio refers to the ratio between 
spherical and elliptical Mj 





M T (r 20 o)/ X 10 i4 M© 






M T (r 50 o/ x 10 ij M o ) 






Cluster JName 


Spherical 


xilliptical 


Ratio 


Spherical 


lilliptical 


Ratio 


A586 


7.3 ± 3.0 


/>J -3.1 


0.97 ± 0.59 


3.0 ± 1.3 


3.1 ± 1.4 


0.97 ± 0.65 


A611 


5.7 ± 1.1 


5.8 ± 1.2 


0.98 ± 0.30 


2.9 ± 0.7 


2.9 ± 0.7 


1.00 ±0.34 


A AO 1 
AOZl 


f: c+2.4 
6 ' 8 -2.5 


/•2_ 24 


C\ Q/1 _i_ r\ 


Z.U ± L.j 


9 9 + 1-3 

2.2_ 14 


°- 91 -i.oo 


A / / J 


J.l ± 1. / 


/.4 ± Z.j 


n aq -i- n (sf* 
u.oy ± u.oo 


z.4 ± u.y 


ri i +1.4 
3 ' 1 -1.3 


a 77+0.76 

0.77_ 073 


A 7Q 1 

A/ol 


C Q , 1 1 

j.y +1.1 


/.Z ± l.o 


U.oZ ± U.jo 


o o _i_ n a 
Z.y ± U.o 


3.Z ± l.U 


n qi _i_ n a 1 

u.vi ± u.4i 


A990 


9 Q+0.6 

^• y -o.i 


2.9 ± 0.6 


1 oo +0 - 29 

i - -0.21 


1.6 ±0.3 


1.6 ±0.3 


1.00 ±0.27 


A1413 


5.7 ± 1.4 


5.8 ± 1.5 


0.98 ± 0.36 


2.7 ± 0.9 


2.8 ±0.8 


0.96 ± 0.46 


A1423 


3.1 ± 1.1 


4.3 ± 1.8 


0.72 ± 0.76 


1.6 ±0.6 


2.0 ± 0.9 


0.80 ± 0.73 


A1758a 


5 9 +1 -° 


6.2 ± 1.2 


95 +0 - 27 

U ^ J -0.28 


3.6 ±0.6 


3.7 ±0.7 


0.97 ± 0.26 


A1758b 


6.3 ± 2.7 


5.8 ±2.5 


1.09 ±0.56 


3.1 ± 1.4 


2.8 ± 1.2 


1.11 ±0.56 


A2009 


6.6 ±2.1 


5 7 +2 - 6 

J -2.9 


1 16 +0 - 48 

i,iO -0.52 


9 Q+0.3 


2 +L2 


1 45 +0 - 42 

^ -0.53 


A2111 


6.0 ± 1.3 


6.0 ± 1.4 


1.00 ±0.32 


2.6 ± 0.7 


2.7 ± 0.9 


0.96 ± 0.44 


A2146 


7.1 ± 1.0 


7.5 ± 1.5 


0.95 ± 0.26 


3.9 ±0.7 


3.8 ±0.8 


1.03 ±0.27 


A2218 


8.7 ± 1.3 


9.0+|-6 


97+ - 24 


3.9 ±0.9 


4 +1 ° 


97+0.35 

uy 1 -0.33 


RXJ0142+2131 


5 3 +L6 


5 4 +i:s 

3 -1.9 


98 +0 - 46 

U ^°-0.49 


2.4 ± 1.0 


2.5 ± 1.0 


0.96 ± 0.6 


RXJ1720+2638 


2.9 ± 0.6 


3.6 ±0.7 


0.81 ±0.35 


1.7 ±0.3 


2.1 ±0.4 


0.81 ±0.32 



2000). We acknowledge the use of NASA's Sky View facility 
(http://skyview.gsfc.nasa.gov) located at NASA Goddard Space 
Flight Center. 



REFERENCES 

Allen S. W., Fabian A. C, 1998, MNRAS, 297, L63 
Allen S. W., 2000, MNRAS, 315, 269 

Allen S. W., Evrard A. E., Mantz A. B., 2011, ARA&A, 49, 409 
Allison J. R., Taylor A. C, Jones M. E., Rawlings S., Kay S. T., 

2011, MNRAS, 410, 341 
AMI Consortium, et al., 2011, MNRAS, 415, 2708 
AMI Consortium, et al., 2011, MNRAS, 415, 2699 
AMI Consortium, et al., 2012, MNRAS, 419, 2921 
AMI Consortium: Olamaie M. et al., 2012, MNRAS, 421, 1136 
AMI Consortium: Rodriguez-Gonzalvez C et al., 2011, MNRAS, 

414, 3751 ^ 

AMI Consortium: Shimwell. T. W., et al., 2010. larXiv: 1012.4441 1 
AMI Consortium: Zwart J. T. L. et al., 2008, MNRAS, 391, 1545 
AMI Consortium: Zwart J. T. L., et al., 2011, MNRAS, 418, 2754 
Arnaud M., Pratt G. W., Piffaretti R., Bohringer H., Croston J. H., 

Pointecouteau E., 2010, A&A, 517, A92 
Baars J., Davies R., Jorgensen I., Bergmann M, Crampton D., AJ, 

130, 445 

Barrena R., Boschin W., Girardi M., Spolaor M., 2007, A&A, 467, 
37 

Bauer, F, E., et al., 2005, MNRAS, 359, 1481 
Bautz M. W., et al., 2009, PASJ, 61, 1117 
Birkinshaw M., 1999, Physics Reports, 310, 97 
Birkinshaw M., Gull S. F. & Northover K., 1981, MNRAS, 197, 
571 

Birkinshaw M., Gull S. F, 1984, MNRAS, 206, 359 

Birkinshaw M., Hughes J. P., 1994, ApJ, 420, 33 

Bohringer H., Schuecker P., Guzzo L. et al., 2004, A&A, 469, 363 

Bonamente M., et al., 2005, ApJ, 614, 56-63 

Bonamente M., Joy M. K., LaRoque S. J., Carlstrom J. E., Reese 

E. D., Dawson K. S., 2006, ApJ, 647, 25 
Bonamente M., et al., 2012, NJPh, 14, 025010 
Browne I. W. A., Wilkinson P. N., Patnaik A. R., Wrobel J. M., 

1998, MNRAS, 293, 257 



Browne I. W, Mao S., Wilkinson P. N., Kus A. J., Marecki A., 

Birkinshaw M., 2000, SPIE, 4015, 299 
Buote D. A., 2001, ApJ, 553, L15 
Carlstrom J. E., Joy M., Grego L., 1996, ApJ, 456, 75-78 
Carlstrom J. E., Holder G. P., & Reese E. D. 2002, ARA&A, 40, 

643 

Carlstrom J. E., et al., 2011, PASP, 123, 568 
Cavagnolo K. W, Donahue M., Voit G. M. et al., 2009, ApJS, 182, 
12 

Cavaliere A. & Fusco-Fermiano R., 1978, A&A, 70, 677 
Corless V. L., King L. J., Clowe D., 2009, MNRAS, 393, 1235 
Cypriano E. S., Sodre Jr., Kneib L. et al, 2005, ApJ, 630, 38-49 
Donnarumma A., et al., 2010. larXiv: 1002. 1625 1 
David L. P., Kempner J., 2004, ApJ, 613, 831 
Ebeling H., Edge A. C, Bohringer H. et al., 1998, MNRAS, 301, 
881 

Ebeling H., Edge A. C, Allen S. W. et al., 2000, MNRAS, 318, 333 
Eke V. R., Cole S., Frenk C. S., 1996, MNRAS, 282, 263 
Evrard A. E., 1997, MNRAS, 292, 289 

Fabjan D., Borgani S., Rasia E., Bonafede A., Dolag K., Murante 

G. , Tornatore L., 2011, MNRAS, 416, 801 
Feroz F, Hobson M. P., 2008, MNRAS, 384, 449 

Feroz F, Hobson M. P., Bridges M., 2009, MNRAS, 398, 1601 
Feroz E, Hobson M. P., Zwart T. L. et al., 2009, MNRAS, 398, 
2049 

George M. R., Fabian A. C, Sanders J. S., Young A. J., Russell 

H. R., 2009, MNRAS, 395, 657 

Giovannini G., Tordi M., Feretti L., 1999, NewA, 4, 141 

Girardi M., Mezzetti M., 2001, ApJ, 548, 79 

Govoni F, Feretti L., Giovannini G., Bohringer H., Reiprich T. H., 
Murgia M., 2001, A&A, 376, 803 

Govoni F, Markevitch M., Vikhlinin A., VanSpeybroeck L., Fer- 
etti L., Giovannini G., 2004, ApJ, 605, 695 

Govoni F, Murgia M., Markevitch M., Feretti L., Giovannini G., 
Taylor G. B., Carretti E., 2009, A&A, 499, 371 

Grainge K., Jones M., Pooley G., Saunders R., Edge A., 1993, MN- 
RAS, 265, L57 

Grainge K.G.B., et al., 1996, MNRAS, 228, 17 

Grainger W.F., et al., 2002, MNRAS, 337, 1207 

Haines C. P., Smith G. P., Egami E., Okabe N., Takada M., Ellis R. 



© 0000 RAS, MNRAS 000, 000-000 



44 AMI Consortium: Rodriguez-Gonzdlvez et al. 



S., Moran S. M., Umetsu K., 2009, MNRAS, 396, 1297 
Hanisch R. J., 1980, AJ, 85, 1565 

Henriksen M, Wnag Q. D. & Ulmer M.R, 1998, AJ, 116, 1529 
Herbig T., Lawrence C. R., Readhead A. C. S., Gulkis S., 1995, 

ApJ, 449, L5 
Hincks A. D., et al., 2010, ApJS, 191, 423 
Hoshino A., et al, 2010, PASJ, 62, 371 

Hudson D. S., Mittal R., Reiprich T. H., Nulsen R E. J., Andernach 

H., Sarazin C. L., 2010, A&A, 513, A37 
Jones M. E, Saunders R., Alexander R et al., 1993, Nature, 365, 

320 

Jones M. E. et al., 2005, MNRAS, 357, 518 
Joy M., et al., 2001, ApJ, 551, LI 
Kaiser N., 1986, MNRAS, 222, 323 

Katayama H., Hayashida K., Takahara E, Fujita Y, 2003, ApJ, 585, 
687 

Kawaharada M., et al., 2010, ApJ, 714, 423 

Klein U., Rephaeli Y., Schlickeiser R. et al., 1991, A & A, 244, 43 

Korngut R, Dicker S., Reese E. D., Mason B. S., Devlin M. J., 

Mroczkowski T., Sarazin C. L., Sun M., 201 1, A AS, 43, #227.05 
Korngut P. M., et al., 2011, ApJ, 734, 10 
Kotov O., Vikhlinin A., 2005, ApJ, 633, 781 
Krause E., Pierpaoli E., Dolag K., Borgani S., 2012, MNRAS, 419, 

1766 

Kravtsov A. V., Vikhlinin A., Nagai D., 2006, ApJ, 650, 128 
Lancaster K., Birkinshaw M, Gawronski, Browne I, Feiler R. et al., 

2007, MNRAS, 378, 673 
Lancaster K., et al., 2011, MNRAS, 418, 1441 
LaRoque S. A., Bonamente M., Carlstrom J. E., Marshall K. J. et 

al., 2006, ApJ, 652, 917-936 
Larson D., Dunkley J., Hinshaw G., Komatsu E. et al., 201 1, ApJS, 

192, 16 

Lo K., Martin R., Chiueh T., 2001, aprs.conf, 235 
Partridge R. B., Perley R. A., Mandolezi N., 1987, ApJ, 317, 112 
Pratt G. W., Bohringer H. & Finoguenov A, 2004, A&A, 429, 
791-806 

Machacek M. E., Bautz M. W., Canizares C. et al., 2002, ApJ, 567, 
188 

Mahdavi A., Hoekstra H., Babul A., Henry J. P., 2008, MNRAS, 
384, 1567 

Mandelbaum R., Seljak U., 2007, JCAP, 6, 24 
Markevitch M., 1997, ApJ, 483, L17 

Markevitch M., Gonzalez A. H., David L., Vikhlinin A., Murray S., 

Forman W., Jones C. & Tucker W., 2002, ApJ, 567, L27 
Marriage T. A., et al., 2011, ApJ, 737, 61 
Marrone D. P., et al., 2009, ApJ, 701, LI 14 
Marrone D. P., et al., 2011, arXiv JarXiv: 1107.5 115 1 
Marshall P. J., Hobson M. P. & Slozar A., 2003, MNRAS, 346, 489 
Massardi M., Ekers R. D., Ellis S. C, Maughan B., 2010, ApJ, 718, 
L23 

Maughan B. J., Jones C., Forman W., Van Speybroeck L., 2008, 
ApJS 

Mazzotta P., Markevitch M., Vikhlinin A., Forman W. R., David 

L. P., van Speybroeck L., 2001, ApJ, 555, 205 
Mazzotta P., Rasia E., Moscardini L., Tormen G., 2004, MNRAS, 

354, 10 

Mazzotta P., Giacintucci S., 2008, ApJ, 675, L9 
McGlynn T. A., Fabian A. C, 1984, MNRAS, 208, 709 
Menanteau F, et al., 2011, arXiv JarXiv: 1109.0953 1 
Meneghetti M., Rasia E., Merten J., Bellagamba F, Ettori S., Maz- 
zotta P., Dolag K., Marri S., 2010, A&A, 514, A93 
Miller C. J., et al., 2005, AJ, 130, 968 



Moln ar S. M., H earn N. C, Stadel J. G., 2012, arXiv, 

larXiv: 120 L15331 
Morandi, A., Ettori, S., Moscardini, L. 2007, MNRAS, 379, 518 
Morandi, A., Ettori, S., 2007, MNRAS, 708 
Motl P. M., Hallman E. J., Burns J. O., Norman M. L., 2005, ApJ, 

623, L63 

Mroczkowski T., et al., 2009, ApJ, 694, 1034 

Mroczkowski T., et al., 2011, MmSAI, 82, 485 

Mroczkowski T., 2011, ApJ, 728, L35 

Muchovej S., et al., 2011, ApJ, 732, 28 

Nagai D., Vikhlin A & Kravtsov A., 2007, ApJ, 655, 98 

Nagai D., Lau E. T., 2011, ApJ, 731, L10 

Navarro J. F, Frenk C. S., White S. D. M., 1996, ApJ, 462, 563 

Nelson K., Rudd D. H., Shaw L., Nagai D., 2011, arXiv, 

larXiv: 11 12.36591 
Okabe N., Umetsu K., 2008, PASJ, 60, 345 
Okabe N., Takada M., Umetsu K., Futamase T., Smith G. P., 2010, 

PASJ, 62,811 

Olamaie M., Hobson M. P., Grainge K. J. B., 2011, arXiv, 

larXiv: 1109.28341 
Owers M. S., Nulsen P. E. J., Couch W. J., 2011, ApJ, 741, 122 
Patnaik A. R., Browne I. W. A., Wilkinson P. N., Wrobel J. M., 

1992, MNRAS, 254, 655 
Planck Collaboration, et al., 2011, A&A, 536, A8 
Planck Collaboration, et al., 2011, A&A, 536, A10 
Planck and AMI collaborations et al., 2012, in preparation 
Poole G. B., Babul A., McCarthy I. G., Fardal M. A., Bildfell C. J., 

Quinn T., Mahdavi A., 2007, MNRAS, 380, 437 
Pratt G.W, Arnaud M., 2005, A&A, 429, 791-806 
Lau E. T., Kravtsov A. V, Nagai D., 2009, ApJ, 705, 1129 
Rasia E., Tormen G., Moscardini L., 2004, MNRAS, 351, 237 
Rasia E., et al., 2012, arXiv JarXiv: 1201. 15691 
Richard J., et al., 2010, MNRAS, 404, 325 
Ricker P. M., Sarazin C. L., 2001, ApJ, 561, 621 
Rines K., Geller M. J., Diaferio A., 2010, The Astrophysical Jour- 
nal Letters, 715, 180 
Rizza E., Burns J. O., Ledlow M. J. et al., 1998, MNRAS, 301, 328 
Romano A., et al., 2010, A&A, 514, A88 

Rossetti M., Eckert D., Cavalleri B. M., Molendi S., Gastaldello E, 

Ghizzardi S., 2011, A&A, 532, A123 
Rozo E., Wu H.-Y, Schmidt F, 2011, ApJ, 735, 118 
Rudnick L., Lemmerman J. A., 2009, ApJ, 697, 1341 
Russell, H.R. and Sanders, J. S. and Fabian, A. C. and Baum, S. 

A. and Donahue, M. and Edge, A. C. and McNamara, B. R. and 

O'Dea, C. P., 2010, MNRAS 
Sanderson, A. J. R., Edge, A. C, Smith, G. P., 2009, MNRAS, 398, 

1698 

Saunders R., et al., 2003, MNRAS, 341, 937 

Schmidt R. W, Allen S. W, 2007, MNRAS, 379, 209 

Sehgal N., Hughes J. P., Wittman D., Margoniner V, Tyson J. A., 

Gee P., dell' Antonio I., 2008, ApJ, 673, 163 
Sifon C, et al., 2012, arXiv JarXiv: 1201.0991 1 
Sijacki D., Springel V, Di Matteo T., Hernquist L., 2007, MNRAS, 

380, 877 

Squires G., Kaiser N., Babul A, Fahlman G, Woods D, Neu- 
mann D. M., Bohringer H., 1996, ApJ, 461, 572 

Smith G. P., Edge A. C, Eke V. R., Nichol R. C, Smail I., 
Kneib J. P., 2003, ApJ, 590, L79 

Smith G. P., Kneib J. P., Smail I, Mazzotta P., Ebeling H., 2005, 
MNRAS, 359, 417 

Snowden, S. L., Mushotzky, R. F, Kuntz, K. D., Davis, D. S., 2008, 
A&A, 478, 615 



© 0000 RAS, MNRAS 000, 000-000 



SZ study of 19 LoCuSS galaxy clusters with AMI out to the virial radius 

Sunyaev R. A., & Zel'dovich Y. B., 1970, Comments on Astro- 
physics and Space, 2, 66 

Takizawa M., Nagino R., Matsushita K., 2010, arXiv, 
larXiv: 100433221 

Tauber J. A., et al., 2010, A&A, 520, Al 

Tsuboi M., Miyazaki A., Kasuga T. et al., 2004, PASJ, 56, 711 

Venture, T., et al, 2008, A&A 484, 327-340 

Venturi et al, 2008, A&A, 484, 327 

Vikhlinin A., et al., 2005, ApJ, 628, 655-672 

Voit G. M.,2005, RvMP, 77, 207 

Wang Q. D, Ulmer M. P. & LAvery R. J., 1997, MNRAS, 288, 702 
White S. D. M., Navarro J. R, Evrard A. E., Frenk C. S., 1993, 

Natur, 366, 429 
White D. A., 2000, MNRAS, 312, 663 

Wilkinson P. N., Browne I. W. A., Patnaik A. R., Wrobel J. M., 

SorathiaB., 1998, MNRAS, 300, 790 
Williamson R., et al., 2011, ApJ, 738, 139 
Wu J.-H. P., et al., 2008, MPLA, 23, 1675 

Vikhlinin A., Kravtsov A., Forman W, Jones C, Markevitch M., 

Murray S. S., Van Speybroeck L., 2006, ApJ, 640, 691 
Voit G. M., Kay S. T., Bryan G. L., 2005, MNRAS, 364, 909 
Zhang Y. Y, Okabe N., Finoguenov A. et al., 2008, A&A, 482, 

451 

Zhang Y. Y, Okabe N., Finoguenov A. et al., 2010, ApJ, 711, 1033- 
1043 



© 0000 RAS, MNRAS 000, 000-000 



